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ADVANCE RESTRICTED REPORT 

AN INVESTIGATION OF AIRCRAFT HEATERS 

XVIII - A DESIGN MANUAL FOR EXHAUST GAS 

.AND AIR HEAT EXCHANGERS 

By L. M. K. 3oelter, R. C. Martinelli, 
F. E. Romie, and E. H. Morrin 

SUMMARY 


Heat exchangers for the transfer of heat from one fluid 
to another fluid at, -a lower temperature are "basic elements 
of modern airplane installations. Examples of such ex- 
changers are intercoolers, oilcoolers, cahin-air heaters, 
exhaust gas heat exchangers for wing ant i - i c ing 'sy s t emg , and 
the heated wing itself. The "basic elements for the design 
of oilcoolers and intercoolers are covered by the report en- 
titled "Design, Selection, and Installation of Aircraft Heat 
Exchangers" by George P. Wood and Maurice J. Brevoort. The 
following report is concerned with the elements of design of 
cabin-air heaters, ving anti-icing heat exchangers, and other 
exchangers in wh’ich the hot fluid consists of air (or the 
products of combustion of air and a hydrocarbon fuel), and 
the cold fluid is -air. 

This report, which summarises a series of reports issued 
by the NACA under the title "An Investigation of Aircraft 
Heaters" (I to XXIII), is divided into four parts. The basic 
equations for the determination of the thermal resistances 
involved in heat exchanger design are presented in f part I. 
Several examples of the application of these basic equations 
of part I to the prediction of the thermal performance of a 
number of heaters are. presented in part II. Nonisothermal 
pressure drop, is* discussed in part III, However, isothermal 
pressure drop characteristics are not presented in detail 
since several readily available references cover this aspect 
of heater design. Also in part III the heat requirements of 
aircraft are discussed briefly; the equations used to correct 
heater performance to any altitude are presented; and the 
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equations required to predict the performance of a ram- 
operated heater-and-duct system at any airplane speed and 
altitude are summarized. 

Part IV consists of an appendix in which the physical 
properties of air are given for a range of temperatures from 
-100° to 1600° F, 

INTRODUCTION 


The design of an exchanger to transfer heat from the 
products of combustion of a hydrocarbon fuel to air, with 
the two fluids separated by a metallic wall, presents a com- 
plex problem. A complete design of a heat exchanger system 
must include consideration of. the following data: 


1 . Heat requirements of the system to which the hot air 
is being supplied . The air may go to a cabin, to a wing anti- 
icing system, to the carburetor air intake, and so forth. 

Such items as heat losses through cabin walls, air leakage, 
heat losses along the airfoil, and so forth, must be known 

in order to establish the heater output necessary to perform 
the desired task. * 

2 , The allowable pressure drop through the heater. 

This specification depends upon the rates of air flow desired 
through the heater, the available total-pressure difference 
across the duct system, the design of the duct work, and so 
forth. As experience is gained with heater installation, 
maximum allowable pressure drops can be specified. • Care must 
be taken, in particular, to design the duct work leading to 
and from the heater as carefully as possible. For example, 
in many installations' the bend leading the hot air from the 
heater is poorly designed because of space limitations so 
that the pressure loss across this element of the duct sys- 
tem may be several times the total pressure drop across the 
heater itself. 


3. The heater design may be established once the heater 
output and allowable pressure drops are specified. This re- 
port, in particular, includes the consideration of the ther>- 
mal aspects of the design. The pressure drop discussion is 


♦Methods of making these calculations or references 
thereto are presented in this report. 
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not detailed because these data may be found in various 
readily available publications. (See references 1, 2, 3, 4, 
5,. 6, 7, 8, 9, and 10.) Certain pressure loss phenomena are, 
however, discussed in part III of the report. 

In addition to the thermal and flow analyses of the 
heater such items as available space, heater weight, heater 
life, and so forth, must be given careful consideration be- 
fore the final heater design is established. 

9- 1 

4, Design of Ducting . The design of the duc.t work 
leading the hot air from the heater to the desired location 
should be considered simultaneously with the- heater design. 
The advantage of a low pressure drop heater may be wholly 
invalidated by poorly designed duct work. As mentioned pre- 
viously, isothermal pressure losses in duct work can be es- 
timated from data available in the literature, A discussion 
of the noniaothermal performance of the ducting will be 
found in part III of this report. 

Many of the equations presented in this report require- 
further experimental verification. Research of this nature 
is being carried out now in the Mechanical Engineering Lab- 
oratories of the University of California. 

I 

It is suggested that, before using the . inf ormat i on in 
this report, part I section A be read carefully, since basic 
concepts and definitions are outlined in this part. 

This investigation, conducted at the Mechanical Engi- 
neering Laboratories of the University of California, was 
sponsored by and conducted with the financial assistance of 
the National Advisory Committee for Aeronautics. 


I. SUMMARY OP HEAT TRANSFER EQUATIONS 
A. GENERAL DISCUSSION OP HEAT TRANSFER MECHANISMS 


Heat is transferred by three mechanisms? conduction, 
convection, and radiation.* 


*Heat also may be transported by the process of evapora- 
tion or condensation, but this mechanism is usually termed 
"mass transfer," rather than "heat transfer." 
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Conduct! on may be defined as heat transfer through a 
body, unaccompanied by any appreciable motions of matter. 

Convection may be of two types: free convection and 
forced convection. 

1 Free convection is the phenomenon of heat transfer 
from a stationary body to a "stationary" fluid - 
that is,' a fluid which is at rest except for the 
convection currents set up by the buoyant forces 
resulting from the heating (or cooling) of the fluid 
in immediate contact with the surface of the body, 

In the normal gas-air heater design, free convection 
is unimportant. 

2. Forced convection may be defined as heat transfer 

from a body to a fluid which has a velocity relative 
to the body, the flow being set up by some external 
agency, such as a pump, a fan, or ram pressure. 

Radi at i on concerns energy which travels in the form of 
electromagnetic waves. Light, radio waves, X-rays, radiant 
heat, and so forth, all are forms of radiation. Heat is 
transferred by radiation when radiant energy starts from a 
body excited thermally, travels across an intervening space, 
and is finally absorbed by another body. . 

In the heat exchangers to be discussed in this report, 
in which heat from the product's of combustion of a fuel is 
transferred to air, all three mechanisms may occur simul- 
taneously. 

In. figure 1 is shown a typical temperature distribution 
in the fluid streams of a gas-air heat exchanger. As the 
hot gas passes the heat exchanger surface, heat is transferred 
to it by forced convection. In addition, some of the gases 
in the products of combustion will transfer heat to the sur- 
face by the process of gaseous radiation. * The sura of the 
heat transferred by convection and radiation then must pass 
through the heat exchanger wall by conduction. Finally, the 
heat is transferred by the process of forced convection to 
the. cold air moving along the wall. Some heat may be lost 
also by radiation from the heat transfer surface if it "sees" 
another surfac e at a lower temperature. 

♦The gaseous radiant heat transfer usually is less than 
10 percent of the convective heat transfer and may be neg- 
le c ted in a practical design. (See reference 11.) 
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Temperature, °F 


Figure 1.- Typical temperature distribution in the fluid 
streams of a gas-air heat exchanger. 


1. Conduct i on . - The elementary equation for the unidi- 
rectional conduction of heat in the steady state through a 
solid, sometimes called Fourier's law, may he written in its 
simple form as: 

q = -kA — — ( 1 ) 

dx 

where 

q rate of heat transfer hy conduction, Btu/hr 
k thermal conductivity of the solid, Btu/hr ft 8 
t temperature, which is independent of time, °F 
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temperature gradient in the solid, F/ft (dt/dx is 
negative when the temperature decreases with in- 
. creasing x ) 

distance, measured in direction of heat flow, ft 

area, at any x, through which heat is flowing, meas- 
ured in a plane perpendicular to direction of heat 
flow, ft 2 

For the simple, case of heat flow through a plane wall 
of thickness L, with temperature independent of time, and 
with flow in the x-direction only, equation (l) may be inte- 
grated : 



If k is independent of t and x, 

q = (t a - t 2 ) = = £t ( 2 ) 

L L R 

where 

R = L/kA "thermal resistance" offered by a plane wall to 
heat transfer by conduction, °F/(Bfu/hr) 

Equation (l) also may he readily integrated for the 
case in which k varies with temperature in a simple manner, 
and for the case of unidirectional heat flow ty -conduction , 
in cylinders, spheres, and other shapes for which A may he 
readily expressed algebraically as a function of x f (See 
reference 12, p. XIV-27,) In most heat exchanger designs 
equation (2) is sufficiently accurate. 


dt/dx 


x 

A 


a* [< 
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Figure 2.- Temperature distribution in wall of thickness 
L. k uniform; t 1 , tg uniform. 


2 . Convection . - The rate of heat transfer by convection 
from a solid to a fluid is controlled by the conduction of 
heat through the fluid immediately in contact with the solid 
surface. The fluid directly in contact wi.th the solid is at 
rest relative to the solid and thus the Fourier law of con- 
duction- can be written for the flow of heat from the body to 
the fluid: 



where 


rate of heat transfer per unit area by convection, 

-Jill ’ 

ft* hr 

B tu 

thermal conductivity of the fluid, r— 

hr ft* 

ft 
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temperature gradient in the fluid immediately in eon 
tact with the solid, °F/ft; the subscript y=o 
refers to the surface of the solid; the partial 
differential form is used to emphasize that a 
particular point on the surface of the body is 
being considered.* 


The temperature gradient 



is controlled by the 


type of flow existing near the surface of the body and by 
the physical properties of the fluid. 


Although equation (3) is strictly the basic equation 
for convective heat transfer, engineers for many years have 
used the well-known "Newton's law of cooling" to express 
heat transfer by convection. This is usually written as! 

(t- - T ) 

q - to A(tg - T J » — — — — (3a) 


where 

q rate of heat transfer by convection, Btu/hr 

3 

A heat transfer area, ft 




surface temperature of solid, °F 


T 


00 


*c 


R = 


1 

f c A 


temperature of fluid far from the body, °F 

unit thermal convective conductance (sometimes 

called the film transfer factor or heat transfer 
coefficient), Btu/hr ft a F 

thermal resistance to convective heat transfer, be- 
tween the temperature t a and , °F/ (Btu/hr) 


Equation (s) defines the unit thermal conductance f c . 


♦Throughout this report the lower case (roman) letter 
t is used to signify the temperature of the solid, and the 
Greek letter T to signify the temperature of the fluid. 
Capital T signifies absolute temperatures of either the 
fluid or the solid. 
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The unit thermal convective conductance fc is usually 

determined experimentally "by measuring q, A, t 8 , and T eo 
and is correlated for use by designers by means of various 
dimensionless moduli to be mentioned. 

It is, instructive to note that equations (3) and (3a) 
are both expressions for heat transfer by convection. By 
equating the two, a basic definition of fc can be estab- 
lished. Thus: 


Of),..' 

- t j 


It is convenient to express equation (3b) in dimension- 
less form: 


t 



(3c) 


where l is a significant dimension in the system which 
fixes the geometry of the solid object from which convection 
is occurring. This point is discussed further. 


f c 1 ' 

The dimensionless modulus — : is called the Nusselt 

k f 


number in heat transfer work. It is noted that the Nusselt 
number is equal to the temperature gradient in the fluid im 
mediately in contact with the solid, divided by the ratio 



From this strict definition of the Nusselt num- 


l 


her several conclusions can be drawn. 

var i e s over . 

y=o 

the surface of a solid, fg, and thus the Nusselt 
number, will also vary from point to point. Most 
heat transfer experimenters to date have contented 
themselves with measurement of the average Nusselt 
number. In this report an effort is made to 
emphasize the point variation of the Nusselt number. 


1, Since the temperature gradient 
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2, Any sharp change in temperature along the surface of 

the solid \irill produce sharp changes in f c and 
thus the Nusselt number. 

3. The Nusselt number will be a function of the factors 

which determine the temperature gradient in the 
fluid immediately in contact with the 90 lid. 


It 


strated 



has been shown analytically and has also been demon- 
experimental ly that the temperature gradient 

is mainly a function of the flow conditions which 


exist next to the object from which heat is being lost by 
convection. (See reference 13.) It is well known that in 
steady state forced convection the flow conditions are char- 
acterized by the Reynolds number of the flow system and the 
shape (including the roughness) of the object. The Reynolds 
number may be written as! 


He . Mil . (t) L_ . ... 9»._ 

US W 3600 lig 3600 |.g 

where J 


(4) 


l a significant dimension of object over which the fluid is 
flowing, ft 

(The choice of this dimension depends upon the geomet- 
rical system being considered. Thus at the entrance 
to a pipe, the significant dimension is the distance 
from the entrance of the pipe to the point under con- 
sideration; once the entrance section has been trav- 
ersed, however, the significant dimension becomes the 
pipe diameter. It is obvioxis that the significant di- 
mension always must be specified when the Reynolds 
number is stated.) 

M- absolute viscosity of fluid, lb sec/ft 2 
u velocity' of fluid, ft/sec 
y weight density of fluid, lb/ft 3 


(Since in many cases, u, p. 
point in the fluid stream, 
number must be accompanied 
by which the magnitudes of 
lished in the calculations 


and y vary from point to 
the statement of the Reynolds 
by a designation of the manner 
u, p., and Y were estab- 
of Re.) 
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g gravitational force per unit 

W fluid flow rate, lb/hr- 

A cross-sectional area through which fluid is flowing, ft 2 

G weight rate of flow per unit cross-sectional area (w/a), 
lb/hr ft 2 

The Reynolds number is representative of the ratio of 
acceleration forces to viscous forces in the fluid stream 
and therefore is a nondimensi onal parameter. In particular, 
small* magnitudes of the Reynolds number signify viscous 
flow, in which fluid particles flow parallel to each other 
with practically no mixing. Large* magnitudes of the 
Reynolds number indicate turbulent flow, during which ap- 
preciable mixing of the fluid occurs due to the eddies and 
vortices resulting from the instability of the turbulent 
motion. Because of the large magnitudes of the Reynolds 
number utilized in normal heater operation, .except for the 
viscous flow near the leading edges of airfoil sections, 
cylinders, and flat plates, turbulent flow v/ill be assumed 
to exist in all heater designs described in the remainder of 
the report. 


of mass , 32.21b 


/(^) 


' the 


In addition to the Reynolds 
system, which establish the 


number and the geometry of 
flow pattern (velocity dis- 


tributi on ) , 


the temperature gradient 



and thus the 


Nusselt number is also 
the fluid, such as its 
specific heat. It has 


a function of certain properties of 
thermal conductivity, viscosity, and 
been shown both by theory and experi- 


*The terms "small" and "large" should be regarded in a 
relative sense, since the numerical magnitude of the critical 
Reynolds number associated with changes from viscous to tur- 
bulent motion depends upon the geometrical arrangement over 
which the flow is taking place, A nominal magnitude for the 
critical Re (based on pipe diameter D) for flow inside 
pipes is 2000, For flow across cylinders, however, the mo- 
tion on the front side of the cylinder may be laminar for 
Re (based on pipe diameter) as high as 10 5 . The flov; in 
the wake behind the cylinder, on the other liana, may become 
turbulent for Re as low as 1800. (See reference 14, pp. 

418 and 423.) (Viscous eddies may form behind a cylinder 
for Re as low as 1.0.) 
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ment that these fluid properties enter the problem as a di- 
mensionless group called the Prandtl number. The Prandtl 
number can be expressed as: 


Pr 

M- 

C P 

g 

k 


3600 |i c p g 

■ Prandtl number (dimensionless) 

k 

/ 2 

absolute viscosity of fluid, lb sec/ft 

o 

heat capacity of fluid Btu/lb F 
gravitational force per unit mass, 

32 - 2i V(“) - 

thermal conductivity of the fluid, 
,Btu/hr 


Thus at each point along the surface of the solid there may 
be written: 



where the subscript x refers to the point values of Nu 
and Re. 


Usually, in the literature, the average Nu is given 
rather than the point values, thus 


Nu 


av 



( 5a ) 


In numerous cases, for many experimental data, equations 
(5) and (5a) can be expressed as a simple power function, 
which will apply in a limited range of the variables. For 
example, for turbulent flow. in smooth, long pipes (the ef- 
fect of the variation of f c at entrance being negligible), 
it has been found that the following relation allows the pre- 
diction of the average fc to be made successfull y: * 

''This equation, discussed in reference 15, is for all 
practical purposes equivalent to that presented by McAdams. 
(See reference 16, p, 168 of 2d ed. ) 
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f c D 


= 0.022 


( 


GD 

36~00” p,g 


)-( 


3600 |i c p g 


) 


0.333 


( 6 ) 


The form of equation (6) is very useful in that it apt- 
plies to any fluid flowing through long tubes in turbulent 
motion. Equation (6), however, is unnecessarily clumsy to 
handle when a particular fluid is under consideration. Since 
this report is concerned mainly with heat exchangers in which 
air is the fluid, equation (6) may be greatly simplified by 
expressing each property of air* ** appearing in the equation 
as a function of temperature, and combining the resulting 
fimctions into a power function of the absolute temperature. 
(See reference 15.) This method has been followed throughout 
this report, so that, although all the equations presented 
are based on generalized forms such as those expressed by 
equation (5), the final form involves only the absolute tem- 
perature, significant dimensions, and the weight rates per 
unit area. For example, by the application of the method 
just discussed, equation (6) for air is reduced to the for 


f c = 5.4 x 10 


-4 


T °' 3 Q°’ 8 

T,0. 3 


(7) 


*Jhe properties of air also may be used for exhaust 
gases with fair accuracy. Reference 17 shows a very small 
difference between the viscosity of exhaust gases and air. 
The thermal conductivity of exhaust gases* is not well known 
at high temperatures. The heat capacity of exhaust gases 
may be calculated if the composition of the gas is known. 
(See reference 18.) 


** In earlier reports of this series the exponent of T 
in several equations was given as 0,296. In order to avoid 
giving the impression of undue accuracy, this exponent has 
been changed to 0.3. For each equation, then, the coeffi- 
cient of T has been changed also in order to give the same 
numerical value as before. In other words, the exponent of 
T has been slightly increased and the coefficient has been 
slightly decreased, so that the result has not changed. 
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where 

f c unit thermal convective conductance 
T mixed-mean absolute temperature of fluid, °R 
G weight' rate per unit cross-sectional area, lb/ft 2 hr 
D tube diameter, ft 

This form of equation is very easy to use and yields results' 
which agree with the more general form (equation (6)) within 
about 2 percent. 

3. Ra d i at i on . - When a body is heated, it emits radiant 
energy at a rate dependent upon its absolute temperature. A 
Planckian radiator - that is, a body that absorbs all the 
radiant energy incident upon it (sometimes called a black 
body) - emits radiant energy at a rate proportional to the 
fourth power of its absolute temperature T a , The rate of 
radiant energy transfer from a Planckian radiator radiating 
to evacuated space at absolute zero is (reference 16, ch. 

Ill of 2d ed. ) 



where 

q r rate of'heat transfer by radiation, Btu/hr 

o S t e f an-B ol t zman radiation constant 

0.173 Btu/hr ( ft 2 
VLOO/ 

A r area ox" body,* ft 2 

If the Planckian radiator, instead of radiating to space 
at absolute zero, radiates to surroundings (also Planckian, 
at a temperature T s ) which completely surround the body, 
the net rate of radiant energy transfer is given by 


*The area of the radiating body effective in radiation 
may not be equal to the total surface area of the body if 
the surface of the body has "re-entrant angles* 11 (See pt. I, 
sec* G- for a further discussion of this point.) 
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= a A 


r r^Y . plY l 

r L Vioo/ \iooy J 


O) 


No actual substance meets the specifications- of a 
Planckian rad. ator, but as a first approximation (reference 
16, ch'. Ill of 2d ed. ) some substances may be considered 
"gray" - that is, they may be defined as absorbing the same 
fraction of the radiant energy incident upon them at all 
wavelengths. -The rate of heat transfer from a small "gray 
body" at temperature Tj to a Planckian body which com- 
pletely surrounds it at temperature T s , is 


ir = [ (ri o) ■ (too) . 


( 10 ) 


where is the emissivity of the radiating gray body 

(always less than unity). If the hody at the temperature T 2 
is not a Planckian radiator and if the two bodies possess a 
given geometrical relationship to each other, the rate of 
heat transfer is given by 


- crA r F AE 


r /t, v n 


O 

i 

fo v 

o 


(IT) 


where F^ is a modulus which modifies the equation for the 

net radiation between Planckian radiators to account for the 
emissivities and relative geometry of the two bodies. (See 
reference 16, pp. 54-60 of 1st ed. ; also references 19 and 
20.) Several values* for F^ are given in part I, section 
G. 


In many engineering applications, a body at temperature 
t lf °F ( T 1 , °R) loses heat by radiation to the surroundings 

at a temperature t 2 , °F (T 2 , °R) and at the same time loses 

heat by convec.ti on to a surrounding gas at a temperature T a . 

In order to simplify this problem, an equivalent unit thermal 
conductance for radiation f r can be defined by an equation 

similar to that used for the definition of unit thermal con- 
ductances for convection. Thus 


*A mechanical integrator (reference 21 ) can be used as 
an aid in determining the shape modulus 
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Figure 3.- Combined radiation and convection. 


q r = f r A (t a - T a ) (12a) 

where A is the total heat transfer surface area of the 
radiating body. After the rate of radiant heat transfer is 
written in this form, the radiant and convective heat trans- 
fer rates may be added as shown; 

q T = q r + q c ~ ^ f r + f c> A “ T a') (12b) 

The use of the equivalent conductance f r reduces the ra- 
diant heat transfer equations to the form of Newton's law 
of cooling.* In order that the radiant heat rate given by 
equation (12a) be equal to that expressed by the more funda- 
mental equation (ll), the equivalent unit conductance for 
radiation f r must be defined as 

*7t should be' noted, that if the body gains heat by con- 
vection and loses heat by radiation (or vice-versa) which is 
often the case, the unit conductances f c and f r in equa- 
tion (12b) will be of oppos i t e sign. In particular, if a body 
is gaining: as .much heat by convection aa it is losing "by ra- 
diation, so that q T = °» the two unit conductances f r and 

f_ are equal and opposite in sign. (See Example, fig. • ) 
c 
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cA r F AE 


= 


(ii.) 4 . (ii_y i 
\i 00/ ViQQ L J 


A (tj - T a ) 


(13) 


Note that T 1 = ( t , + 460) 

The equivalent conductance f r therefore is a function of 

1. The emissivitiea of the radiating surfaces 

2. The relative geometry of the radiating surfaces 

3. The temperatures of the radiating surfaces and of 

the s u r r ounding fluid 


In the preceding equations the following symbols were uti- 
lized: 

a 

A total surface area of heat transfer, ft 

s 

A effective surface area of the radiating body,** ft 
r 

z o 

f_ unit thermal conductance for convection, Btu/hr ft F 
c 

f equivalent unit thermal conductance for radiation, 

Btu /hr ft 3 °F 

F aE combined shape arid emissivity modulus defined by equa- 
tion (8) (sometimes set equal to x Fjj), dimen- 

sionless 

F/^ shape modulus to account for relative geometry of sur- 
faces exchanging heat by radiation (dimensionless) 

Fjj emissivity modiilus to account for emission characteris- 
tics of surfaces ‘ (dimensionless ) 

q heat transfer rate, Btu/hr 

t x temperature of surface . 1, °F 

T\ absolute temperature of surface 1, °R 

*The usual definition of f r does not involve the 
fluid temperature T a . (See reference 16, p. 63 of 2d ed.) 
* f See pt. I, sec, G for a discussion of this term. 
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t 3 
a 

T 

a 


temperature of surface 2, °F 
absolute temperature of surface 2, °R 
St efan-B ol t zmann constant, 0,173 Btu/hr 
temperature of ambient fluid, °E 



ft s 


4, Combined heat tr ansfer mechanisms .- In the previous 
sections the three mechanisms of heat transfer have been con- 
sidered separately. In practice, however, two or more of 
the mechanisms usually occur simultaneously. The transfer 
of heat from one gas to another through a. metallic surface 
illustrates this point. On the hot gas side heat is trans- 
ferred by convection and radiation (reference ll) from the 
hot gas to the wall surface. The total rate of heat trans- 
fer from the hot gas to the wall surface is given by; 


or 







q.T 





Since the steady state exists, the same rate of heat trans- 
fer occurs through the metallic wall by conduction. Thus 

kA / \ 

q T - L l^w g ~ t ^ a ) 

After passing through* the metallic wall, the thermal energy 
is transferred to the cold air hy convection, assuming the 
radiant heat transfer to the air to he negl igihle . * Thus: 


♦Radiant heat transfer to the cold air may he promoted 
hy suspending metal plates in the air stream in view of the 
hot heat exchanger surface. Heat is transferred to these 
plates from the hot surface hy radiation; then heat is trans- 
ferred from the plates to the air stream hy convection. These 
"irradiated convectors" may have an appreciable effect on the 
heat exchanger output hut usually at the expense of greatly 
increased frictional pressure drop. (See example in pt. I, 
sec. G for details.) 
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q.rp - 




By eliminating the intermediate temperatures, t w and t w , 

& & 

the following equation is obtained 



(14) 


The same result could have been readily obtained by 
noting that unidirectional thermal and direct current elec- 
trical circuits are analogous and the corresponding equations 
for the rate of heat flow and current are similar for the 
steady state. (See reference 22.) If the rate of heat flow 
is designated as a thermal current and the temperature dif- 
ferences act as potentials, the following analogous relations 
may be written for the steady unidirectional state. 


AE 


i = (electrical) 


Thus 


.( 


ance s * 


\ kA 

ii- 



f c A. 



q = 


- — — (radiation) 

(~) 



ri) 


may 'be termed "thermal resist- 


The mechanism of heat flow from the hot to the cold gas 
then can he visualized as analogous to the flow of current in 
a simple electrical d irect.-cixrrent circuit. 
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Then, "because 


i = 


Ei «• E4 


i- ♦ A- 

Rj R 2 


+ R, + R, 


then, by analogy, 


T - T 

g a 


The term 


4 


i 


+ L + 

1 





+ f c 
' g 

)\ 

kA 

fc A 



in steady 
presented 

heat 
in th 

flow can he 
e foregoing 

readily analyzed by 
paragraph . 


- 


i 




(15) 



i 

__ -f 

L + _ 

1 



(% + 


A 

kA f 




u 

over-al 1 

conductance UA 

of the thermal 

system 


If the temperature at either surface of the metal is 
desired, inspection of figure 4 reveals, from the thermal 
circuit, that 



( 16 ) 


Then : 



T 


g 


- T 


a 


L t 
kA 
1 

UA ■ 
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or 





) 


Thus the temperature of the metallic surface may be 
estimated readily once the^unit thermal conductances, the 
over-all thermal conductance, and the two gas temperatures 
are known.* It is evident from eqiiation (16) that if the 
ratio (UA/f c A) is small, then the metal temperature will be 
almost equal to the air temperature T a . 


In many cases the values of the thermal conductances 
(f c and U) vary throughout the exchanger. In order to ob- 
tain the temperature of the exchanger surface, the 1 ocal 
values of the thermal conductances at the point in question 
should be used in the foregoing equations. 

The thermal conductances f c may be adjusted to yield 
the proper metal temperatures by changing the fluid velocity, 
and so forth. The effect of these changes, of course, can 
be calculated by means of the foregoing equations. 

In the foregoing equations the symbols have the follow- 
ing significance; 

rate of heat transfer, Btu/hr 

total (convective and radiation) heat transfer rate, 
Btu/hr 

unit thermal conductance for convection on the hot-gas 
side of the heat transfer surface, Btu/hr ft 3 F 

equivalent unit thermal conductance for radiation from 
certain constituents of the hot exhaust gases to the 
heat transfer surface (defined by equation (13)) Btu/ 
hr ft 3 °F. (See reference 11.) 

unit thermal conductance for conveotion on the cold-air 
side of the heat transfer surface, Btu/hr ft 3 °F 

surface area of heat transfer, ft 2 

mean temperature of hot gases, °F 

*In many cases conduction of heat along the metal sur- 
face will have an important effect on the metal temperatures. 
Calculations considering the effect of convection to the metal 
surfaces and also conduction along the metal can be made using 
the "Southwell relaxation method." (See reference 81.) 


q 

q j 



A 

T 
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temperature of heat transfer sur-face in con- 
tact with hot gas, °F 


'W. 


T a 

k 


temperature of heat transfer surface in con- 
tact with cold air, °F 

mixed-mean temperature of cold air, °F 
thermal conductivity of metallic heat trans- 
fer surface, Btu/hr ft 2 



L thickness of heat transfer surface, ft 

E electromotive force, volts 

i current, amperes 

R electrical resistance, ohms 

Ri , R2 , R 3 , R 4 electrical resistances (shown in fig. 4), 

ohms 

E 3 , E 4 voltages (shown in fig. 4), volts 


Mean Temperature Difference 

The rate of heat transfer between the hot and cold gas 
can he written. as; 


1 * DA < T g - T a> 


In the analysis presented in the last paragraph the tem- 
perature difference (t^ - T & ) has been assumed constant. In 

actual exchangers the temperature difference between hot and 
cold gas varies throughout the exchanger. A discussion of 
the proper mean temperature difference to be utilized for 
such cases is presented in part II of this report. 

If the over-all conductance UA and the mean tempera- 
ture difference (t^ - T a} m are known, the thermal performance 
of the heater can be established. 

The remainder of this section presents equations for the 
determination of the unit thermal conductances for the most 
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common flov/ systems met in practice. Combination of the unit 
thermal conductances for any heater in the form of thermal 
resistances, will allow the evaluation of the over-all con- 
ductance UA and thus the heater performance. 


3. FORCED CONVECTION ALONG PLATES 
General 


When a fluid flows along a fixed solid boundary, the 
fluid in immediate contact with the surface attains zero 
velocity. For fluids with relatively low viscosity, such 
as air, the velocity will change from zero to its free stream 
value within a thin layer of fluid next to the wall, called 
the "boundary layer." (Reference 14, p. 50.) The thickness 
of the boundary layer at any point on the body depends upon 
the kinematic viscosity of the fluid, the velocity in the 
free stream, and the geometry of the body. 

If, for example, a flat plate is placed with its sur- 
face parallel to the direction of flow, as the fluid comes 
in contact with the surface of^the plate it is brought to 
rest. At the leading edge of the plate the boundary layer 
is. very thin, since only the fluid in immediate contact with 
the plate has been brought to rest, while the remaining 
fluid flows on with the free-stream velocity u,^. As the 
fluid proceeds along the plate the tangential stresses set 
up by the solid boundary cause more and more of the fluid to 
be retarded, and thus the thickness of the boundary layer 
continually increases. The growth of this boundary layer is 
illustrated in figure 5, in which the velocity distribution 
in the. fluid stream at various points along the plate are 
shown, the vertical scale being greatly enlarged for clarity. 

For this so-called laminar boundary layer, existing near 
the leading edge of the plate, it has been shown that the 
thickness of the boundary layer, the frictional drag on the 
plate (reference 14, p. 50), and the unit thermal conductance 
at any point on the plate (reference 14, p. 623) are all 
functions of the square root of the Reynolds number for the 

( UoPcYN 

plate, l He = ). At a certain point along the surface the 

\ M.g / 

fluid in the boundary layer becomes turbulent, and, from 
this point on (fig. 6), the thickness of the boundary layer, 
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Free-stream velocity* 



Figure 5.- Velocity distribution in laminar 
boundary layer. 
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the drag (reference 1.4, p, 362), and the unit thermal con- 
ductances are functions of the one-fifth power of the 

/ U^xVN 

Reynolds number for the plate (reference 23) (Re = ^ ). 

A number of analytical' solutions for forced convection 
along plates have been presented in the literature (refer- 
ence 14, p. 623), and several investigators have obtained 
experimental data on this flow system (reference 23). On 
the basis of this work, the following equations were derived, 
(See reference 24.) These equations allow the determination 
of the unit thermal conductance for a flat plate along which 
air or exhaust gas is flowing, at any altitude pressure and 
for a temperature range of -60° to 1600° F. 

Laminar Boundary Layer 

(a) Point unit conductance at any point x feet 
from the leading edge of flat plate; 


f c = 0.0562 T f 


0.50 


K x / 


0.50 


(17) 


(b) Average unit conductance in length l (laminar 
sublayer from x = 0 to x = l): 


f„ = T- 


f c . dx 
c x 


= 0.112 T. 


0,50 


(¥) 


0.50 


(18 ) 


Turbulent Boundary Layer 

(a) Point unit conductance at any point x from 
leading edge of flat plate: 

o . e 


f Cv = 0,51 T f 


0.3 (Uoo Y ) 


X 


0 , 2 


(19) 


(b) Average unit conductance in length (turbulent 
boundary layer for x = 0 to x = l): 

l 


■tj 


dx = 0.64 T, 


0 . 3 


(u^Y) 


0 . 8 


l 


0.20 


( 20 ) 
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If a laminar boundary layer exists on part of the plate (fig. 
6) and a turbulent boundary layer over the remainder, the 
average unit conductance in the length may be obtained by 
integration of the point unit conductance. Thus 



/ 


f„ dx + 


f „ dx 


o (laminar ) 


l (turbulent ) 


( 21 ) 


In the foregoing equations, 

f' c point unit conductance at a distance x from leading 
x edge of flat plate, Btu/hr ft 2 °F 

f c average unit conductance of plate in length l, Btu/ 

hr ft 2 °F 

Tf arithmetic average of plate temperature and free air 
stream temperature, °R 

u^ free air stream velocity, ft/sec 

y density of air at temperature Tf , and prevalent pres- 
sure , lb /f t 3 

x distance from leading edge of plate, ft 

distance from leading edge of plate to point at which 
laminar boundary layer becomes turbulent, ft (See 
references 14, p. 361; and 23.) 

I length of plate under consideration, ft 


Example 

The transition from a laminar boundary layer to a tur- 
bulent boundary layer along a flat plate 12 inches in length 
is postulated to occur at a magnitude of the point Reynolds 
number for the plate equal to 50,000. Plate temperature is 
200° F; air temperature, 30° P; free stream air velocity, 

100 feet per second; atmospheric pressure, 14.7 pounds per 

square inch absolute. (a) Calculate the variation of f- 

c x 

along the plate. (b) What is the average f c for the 
plate? 
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(a) Transition point 

The Reynolds number at transition from the laminar 
to the turbulent boundary layer is : 


Re = 


M-g 


= 50,000 


From the data: 


T f = 460 + 2 °° + 30 _ 575 o R 




= 0.406 x 10 " 6 — ■ \ ec (See fig. 42.) 

ft 2 


y = — = 14.7 x 144 _ 0.0692 

RT 53.3 X 575 ft 3 


Thus , 


50000 x 0.406 x 10~ 6 x 32.2 
100 x 0. 0692 


lj - 


li = 0.0943 ft = 1.13 in. 

(b) Point unit conductance 

A laminar boundary layer exists for the first 1.13 
inches of the plate. Thus: 

f 0x = 0. 0562 T f 0 ’ 50 (equation (17.)) 


f c = 0. 0562 X 575 


o,50 A00x0 i 0692 xO * 5 




0 


3,56 B tu 
X V 2 hr ft 3 °F 


On the remainder of the plate there is a turbulent 
boundary layer. Thus for 0.0943 < x < 1.0 feet. 


f c = 0.51 T f 


0.3 (u^V ) 0 * 8 


f Cx = 0-51 (575) 


o .3 (100 x 0.0692) 


- 0.2 

0 . 8 


x °,2 


(equation (19)) 

^- 5 - (0.0943 < x < 1 ft) 

x 0 - 2 hr ft 3 °F 
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C. FORCED CONVECTION INSIDE TUBES 
1. Turbulent Flow in Tubes 


Oeneral . - As air in turbulent motion with a uniform ve- 
locity enters a tube, the fluid immediately adjacent to the 
tube wall is brought to rest. For a short distance along the 
tube wall a so-called laminar boundary layer is formed. (See 
pt. I, sec. B.) If the turbulence in the entering fluid 
stream is sufficiently great, the laminar boundary layer will 
quickly change into a turbulent boundary layer, The turbu- 
lent boundary layer will rapidly increase in thickness until 
it fills the whole of the pipe. From this point on, the ve- 
locity profile across the pipe remains unchanged for isother- 
mal flGw, (See reference 14, p. 360.) 

As the flow pattern changes from a laminar boundary 
layer to a turbulent boundary layer, and finally to the fully 
developed turbulent velocity distribution, the fluid in im- 
mediate contact with the tube wall will maintain its viscous 
motion, altho\igh the fluid further from the wall will flow 
turbulently, (See reference 25.) This viscous layer of 
fluid (initially the laminar boundary layer) is called the 
"laminar sublayer" in the region in which the fully developed 
turbulent velocity distribution has been established. The 
thickness of this laminar sublayer may be readily estimated. 
(See references 26 and 27.) Next to the laminar sublayer a 
"buffer" layer may be visualized in which both viscous and 
turbulent forces play an important part. The remainder of 
the flow system may be considered as a completely turbulent 
"core" in which viscous forces are unimportant, and turbulent 
forces predominate. For a fluid the viscosity of which does 
not change appreciably with temperature, the thicknesses of 
the laminar sublayer and buffer layer are independent of dis^- 
tance along the tube (in the region in which the velocity dis- 
tribution has been established) and depend only upon the 
Reynolds number, based on the tube diameter and the mean ve- 
locity of flow. Figure 8 illustrates these various phenomena. 

♦A very high point unit thermal’ 'conductance exists at 
the entrance to a tube even when a hydr odynami c calming sec- 
tion precedes the heating section, due to the very large tem- 
perature gradient existing at the fluid-wall interface at the 
beginning of the heating section, Latzko in reference 28 
considers this case as well as the case in which the fluid 
enters the tube with uniform velocity. The simplified equa- • 
tions /presented in this section, however, apply more precisely 
to the latter caso. 




Flow Phenomena and Variation in Point C/nil Conductance at the Entrance to a Straight Doc 
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All the equations for the entrance section shown in fig- 
ure 8 are approximations which are 'based, on the postulate 
that the entrance section of the pipe may be considered to 
behave like a flat plate placed parallel to the direction of 
air flow.* (See pt. I, sec. B.) The equations for the unit 
conductance in the portion of the pipe in which the velocity 
distribution is established are close approximations (refer- 
ence 15) to analytical equations which were derived by a 
consideration of the thermal resistance of the laminar sub- 
layer, the buffer layer, and the turb\xlent core (reference 
26). 


In most heater designs, experience has shown that the 
laminar boundary layer occupies a negligible portion of the 
entrance section. If the location of the transition point 
from the laminar to the turbulent boundary layer can be as- 
sumed to exist at the very beginning of the entrance section, 
the laminar boundary layer can be entirely neglected. 

Summary of Equations ; 

Range of Apnlication 

Fluids! air or exhaust gases in smooth tubes 
Temperature range! -60° to 1600° F 
Pressure range! pressure at any altitude 

For short tubes either the point of transition from 
laminar to turbulent boundary layer must be known, or as an 
approximation, a turbulent boundary layer can be postulated 
to exist from the point x = 0, 

*The length of the "entrance section" - i.e., the sec- 
tion at the entrance to the tube in which the velocity and 
temperature distribiit i on have not yet attained their fully 
developed form - has been evaluated in this section by not- 
ing the distance from the tube entrance, at which the point 
unit conductance in the 'entrance section equaled that for the 
remainder of the tube. This procedure (reference 29) yields 
the entrance length l p'4.4 Djj . More precise calculations 

by Latzko (reference 2 b) show that the entrance length is a 
function of the Reynolds number (based on pipe diameter and 
mean velocity). Latzko gives as the entrance length 
l = 0.693 Djj Re 1/4 . At a magnitude of Re = 10,000 this 

equation yields l = 6. S3 Djj , which checks roughly the 
magnitudes of the approximate method. 
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For long tubes the minimum Re = ^6 00 p.g = 

(higher than 2000 to insure fully developed turbulence). 


Short Tubes 


(a) Point unit* conductance (Turbulent Boundary Layer) 
(reference 29) 


f„ = 7.3 x 10 


m 0*3 r <>. 8 

- 4 1 f ^ 


(b) Average unit conductance in length l (Turbulent 
Boundary Layer) 


f c = J / f Cx 4x = 9.1 x 10' 


4 T f 0,3 G°* 8 
l © • 3 


2. Long Tubes 


(& ” 4 ' 4 ) 


(a) Point unit conductance beyond the entrance length 


„ o. e 

f. = 5.4 x 10“ 4 T 0,3 - Q — 

°x PH 


(24) 


(b) Average unit conductance (includes entrance effect) 


in. - , / dx = 5.4 X 10' 


l / * c x 


(25) 


f c point unit thermal conductance, x feet from entrance 
of tube, Btu/hr ft 2 °F 
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f c average unit conductance in length l of tube, Btu/ 
hr fta °F 


x distance from entrance of tube, ft 
l tube length, ft 

Dy hydraulic diameter of tube, equal to 

4 x cross-sectional area 4A 

= — , ft 

wetted perimeter P 

A cross-sectional area of tube, ft 

F wetted perimeter of tube, ft 

G W/A = weight rate of air per unit cross-sectional area, 
lb/hr ft 2 


W weight rate of air, lb/hr 

T arithmetic average of the mixed-mean temperature of the 
• air entering and leaving the tube, °R 

If arithmetic average of air temperature and tube wall tem- 
perature, °R 


Example - Short Duct 

Air at an average temperature of 300° S' flows through 
a tube with the dimensions given in figure 9 at a rate of 
116 pounds per hour. The tube wall temperature is 800° S’, 


,t w = 800 °F 

/ 

W a = 116 lb/hr I 1 

cx 


* a = 300 °F 

4" 



Figure 9.- "Short" duct. 
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wetted perimeter = 4 in. = 0.333 ft 


cross-sectional area = 1 sq in. = 0.00695 ft 2 

hydraulic diameter = — — *_JKJ30695. _ ft 

0.333 


= 0. 0833 ft 


— = 0.333 ft 
12 


l/D H £s_333_ = 4> o (short tube) 
n 0.0833 

Or — 2JJ = 16,700 — 

0,00695 hr ft 2 

_ 800 +_ 300 + 460 _ 1010 ° r 


(a) V/hat is the average unit conductance along the tube? 
% 

(b) What is the variation of unit conductance along the 
tube, assuming a turbulent boundary layer starting at x = 0? 

(a) Check Reynolds number (based on tube hydraulic diameter) 


Re = 


u 16700 X 0. 0833 

3600 ng " 0.50 X 10~ 6 x 32.2 X 3600 


= 24,000 


The viscosity of air at 300° F was obtained from figure 
42. Since Re is quite high, assumption of a turbulent 
boundary layer from x = 0 is probably valid. 

From equation (23) 

f- = 9.1 x 10~ 4 T^° * 3 


f c = 9.1 X 10 


_4 (1010) ( 1 S 7 00 ) ' 

(0.333)°* 2 


21 . 6 


hr ft 3 °F 



Unit conductance 
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Figure 10.- Variation of unit thermal conductance along short duct 
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(t>) Variation of the point unit conductance along the tube. 
From equation (22) 

T 0,3 G 0 * 8 

f = 7.3 x 10 -4 — — = Btu/hr ft 2 °F (x in ft) 

c x x 0 * 2 x° #2 

A plot of f c versus x is shown in figure 10. 

Example - Long Tube 

What is the average unit thermal conductance for the fol 
lowing system (including entrance effect)? 



I length of tube = 3.0 ft 

A cross-sectional area = 4.2 sq in. = 0.0292 ft 2 

P wetted perimeter = 8.4 in. = 0.70 ft 

Dp = -AJL = . 4— X -Q. Q?9? . = 0.167 ft 
n P 0.70 

^ /®H = p ^ 67 = 18 (long tube) 

T average fluid temperature = 200 + 460 = 660° R 
W rate of air flow = 500 lb/hr 

„ 500 . a 

G weight rate per unit area = = 17,100 lb/hr ft 

0 . ©292 

The Reynolds number for the flow is: 



Point unit conductance 


NACA 

ARR 

No. 

5A06 


3! 




Re 

G D 

= 54,800 





3600 u,g 

• 

Thus 

equa 

ti on 

(25) is i 

ippl icabl e 

and , 

0 

f c 

5.4 

o 
1 — \ 

X 

- 4 T 0.3 _G 
D I 

0.8 / 

0.2 ( 1 + 1 
{ v 

•> 4 ) 


R.k X 

icr 4 

x ( 66 o ) 0 ' 3 

.. (17100) 0 * 8 

Yi , 1 1 




A \ LJVJ / 

X CO.167 )°* 3 

ix + x»x ] — _ 

\ 3.0 / hr °F ft 




boundary layer » 


3.0 ft 
8.4 in. 

4.2 sq in, 

17,100 lb/hr ft 2 
200 °F I 



Fully developed 


velocity distribution 


Figure 12.- 


Distance from tube inlet, ft 
Variation of unit conductance along a long tube. 
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A plot of the variation of along the tube length is 

shown in figure 12, The value of fcx the entrance to 

the tuhe was obtained, by means of equation (22), Because 
of. the small effect of the entrance section, magnitudes of 
the mean f c and the point f Cx are practically equal over 

most of the tube length, 

2. Viscous Flow in Tubes 

General . - When a fluid flows through a tube with a 
Reynolds number less than 2000 (based on tube diameter and 
mean velocity), the rate of heat transfer from the fluid to 
the tube walls is very low. The poor heat transfer is a 
direct result of the mechanism of viscous flow,* for in this 
type of fluid motion the particles of fluid tend to follow 
parallel paths with little or no mixing of adjacent layers 
of the fluid. Viscous flow is rarely encountered in aircraft 
heater design, but in some special applications, such as 
heating the interior of the leading edge of a wing by hot 
air, viscous motion of the fluid may result at low rates of 
air flow. 

The determination of the rates of heat transfer during 
viscous motion is complicated by the fact that, due to the 
very slow motion of the fluid through the heat transfer pas- 
sages, appreciable free convection heat transfer may be 
superimposed on the heat transfer by viscous forced convec- 
tion. In references 31 and J2 a method is outlined for the 
approximate determination of the rates of heat transfer due 
to combined free and forced convection in the region of vis- 
cous flow,** when the velocity distribution at the entrance 
to the tube is parabolic and the tube temperature is uniform. 
Based on this work, the following approximate expression wa6 
derived,*** for short vertical circular tubes in which the 
free and forced convection velocities are in the same direc- 

t Ion: 

•See references 16, pp. 189,199 of 2d ed.; 3°» 31# and 
32 for data on viscous forced convection. 

**Bef erence l6, pp. 189,199 of 2d ed. gives a summary of 
average unit thermal conductances in viscous flow in tubes 
and duct 8. 

•••Equation (25a) has been verified experimentally for 
several liquids with a positive coefficient of thermal expan- 
sion flowing vertically upward while being heated, and for 
water being cooled while flowing vertically downward (refer- 
ence 31)* (The viscosity of all the fluids used decreases 
with increased temperature.) 
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4 I 



(25a) 


where 


f c point unit thermal conductance along the tube starting 
from x = 0, Btu/hr ft 3 °F 

D' inner diameter of tube, ft 

V thermal conductivity of the fluid, Btu/hr ft 2 (°F/ft) 

■^W rate of flow of fluid, lh/hr 

c'p heat capacity of fluid, Btu/lb °F 

y 

x distance from entrance of tube, ft 

^ V 

OtT Orashof number (see reference 32 for details) 

/ 

Pr Prandtl number 

A similar equation for heat transfer from two infinite 
vertical parallel plates (uniform in temperature) to a fluid 
being heated while flowing vertically upward (between the 
plates) in viscous motion with a parabolic velocity distribu- 
tion at the entrance to the heating section, is also pre- 
sented in reference 31. Based on this work, the following 
equation may be written, which may be applied to the analysis 
of heat transfer in narrow rectangular ducts (fig. 12a). 



0.98 


59 + 


Wc, 


XX 


8. 

B 


+ 0.153 


^ Gr Pr fiY /4 


(25b) 


where all terms are as defined for equation (25a) except 
6 smallest distance between sides of ducts, ft 
B breadth of duct, ft 
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Simplified Equations 

For use in design of air heat exchangers the following 
further approximations of equations (25a) and (25b) are pre- 
sented.: 

For circular tubes; 




.65 jj- zj 1 + [ 


0.38 W + 3500 D 3 V s At 


(25c) 


For flat ducts : 


f. = 3.80 & 3 i 
°x 5 


1 + 


0.20 W 


3000 § 3 y 3 At 


( 25d ) 


where 

Y average weight density of air flowing through tube, 
lb/ft 3 

At average temperature difference between air and tube walls, 

<>f 


Equation (25c) is applicable for air flowing in round (or 
nearly round) tubes when the Reynolds number for the air 
flow is less than about 2000. Equation (25d) is applicable 
to air flowing viscously in flat ducts. The effect of free 
convection on the rate of heat transfer is approximately ac 
counted for in the equations * ,T the term involving the 
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temperature difference At; the effect of free convection 
usually is small. The equations are most accurate for alti- 
tude pressures and for a temperature near 300° F, hut may he 
used for higher temperatures (1000° F) for approximate values 
of f c . Although equations (25c) and (25d) are based on an 

analysis of vertical tubes and plates with uniform tempera- 
tures and with a parabolic velocity distribution of the 
fluid at the entrance to the heating section, they may be 
used as a rough approximation for tubes and flat ducts with 
•other orientations and with other entrance conditions. The 
equations are least accurate when the velocities aue to free 
convection are in the opposite direction than those due to 
forced convection, 

Exampl e : 

Air flows through a 1-inch. inside diameter round tube 
1.5 feet long under the following conditions.’ 


Mean-air temperature . .. 300° F. 

Air rate 7 lb/hr 

Mean tube surface temperature . . 40° F 

Absolute air pressure . . . . , 10 lb/in 3 


Determine the variation of f c ^ with tube length: 
(a) Calculation of the Reynolds number 


tt 36 00 Dug. 


4 x ? 

3.14 X 3600 ' X X- X 0.498 X 10~ 6 X 32.2 
X o 


■= 1850 


Equation (25c) is applicable, because Re < 2000. 
(b) Fluid Properties (appendix A, Properties of Air) 


Air density V = =■=• = 


RT 53.3 x 760 
Thermal conductivity k = 2.05 x 10“ 2 


10 X 144 _ 0355 113 


ft' 


B'tu 


hr l - 
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Figure 12'b.- Variation of unit con- 



0 .4 .8 1.2 1.6 

Length along duct, ft 


The average f c is obtained by integrating under the 
curve of f c versus x and dividing by l. Thus, for 
round tubes , 
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where 



c = [0.38 w + 3500 D 3 Y s At] 


To simplify the evaluation of f„ for round tubes and 

°av 


ducts, a plot of 


l 



o 


is given in figure 12c as a function of 


l/c. 


Bxample: In the previous example, what is the average f c 

along the tube? 


s o 


Thu s 


f 


c 


av 



3.3 

1.5 ft 

J 0 . 455 = 0.675 

2.05 (fig. 12c) 


3.65 K 
D 


x X 


3, 65 x 0, 0305 
1 

12 


2.05 


1.84 


Btu 
2 0 

hr ft F 


Transiti on Region : 

As has "been discussed, the laminar heat transfer equa- 
tions presented in this section are applicable only below a 
Reynolds number of 2000, and the turbulent heat transfer 
equations are applicable" only ab ove a Reynolds number of 
10,000 # The region between Re = 2000 and Re = 10,000 
cannot be described simply by either the laminar of turbu- 
lent equations, since in this region the fluid flows vis- 
cously in the initial portion of the tube and then breaks 
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Figure 12c.- Function x for use In viscous flow beat transfer equations to obtain average value of f c 
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into turbulent flow in the remainder. Figure 12d represents 
a typical variation of the point unit conductance along a 
tube in which air is flowing with a Reynolds number between 
2000 and 10,000. The point of transition from viscous to 
turbulent flow must be known before the unit conductance in 
the transition region can be determined. 



Figure 12 d.- Typical distribution of f Cx in duct for fluid flowing 
with Reynolds number between 2,000 and 10,000. 


Some evidence exists that at each magnitude of the 
Reynolds number Re (based on tube diameter) transition may 
take place at definite magnitudes of (Re x c/D), but until 
this evidence is investigated further, the evaluation of the 
unit conductance in the transition region can be only a 
rough approximation. 
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D. FORCED CONVECTION ACROSS TUBES 
General 


The flow across a single cylinder can be visualized as 
shown in figure 13. At low magnitudes of the Reynolds number 
(based on the outer diameter of the cylinder and the velocity 
in the free stream) a laminar boundary layer, starting at the 
front stagnation point, exists over the greatest portion of 
the cylinder. Because of the unfavorable pressure gradient 
existing around the cylinder (reference 14, p. 56), separa- 
tion takes place about halfway between the front and rear 
stagnation points, as shown in figure 13. 

The rate of heat transfer, therefore, is controlled by 
a laminar boundary layer on the forward portion of the cyl- 
inder (reference 14, p, 631 ) , and on the rear portion by the 
turbulence existing in the wake of the cylinder The point 
values of f c , measured by Schmidt and Wenner (reference 33) 
on the surface of a right circular cylinder, are shown in 
figure 13. The data reveal that the minimum rate of heat 
transfer exists on the sides of the cylinder. Because of 
the turbulence, at high values of the Reynolds number the 
unit conductance at the rear stagnation point may exceed 
that at the forward stagnation point. 

In the discussion of heat transfer from flat plates 
(pt. I, sec. 3) it was noted that for the laminar boundary 
layer, the unit thermal conductance varied with the 0,5 power 
or the product of velocity and density of the fluid (u^Y) 0,5 , 

For turbulent conditions, the unit conductance varied with 
the 0.8 powet of this product. It is reasonable, therefore, 
that for flow across cylinders the unit thermal conductance 
should vary with a power of (u Y ) between 0.5 and 0.8. Ex- 

OO 

periment ally , for a range of Reynolds numbers (based on tube 
diameter) between 1000 and 50,000, the unit conductance has 
been found to vary with about the 0.6 power of the product. 
(See reference 16, p. 222 of 2d ed. ) 

When several tubes are placed one behind the other, the 
.turbulence produced by the first cylinder will increase the 
heat transfer from the second, and so forth. Thus it will 
be found that the rate of heat transfer from several banks 
of tubes is higher than for the same considered individually/ 
(See reference 8.) 
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A summary of the equations to he used for heat transfer 
for air or exhaust gases flowing across single cylinders and 
tube banks, based on the data of references 8 ; l 4 , p. 63 I; 
l 6 , p, 222 of 2 d ed.; and 33 i® presented. These equations 
une applicable at any altitude pressure and for a tempera- 
ture range of -60° to 1600° P. 

Point Unit Conductances 


i nder 


(a) Unit conductances at front stagnation point of cyl- 


o 4</ u A 0 * 50 

f 0o = °'l 9 4 T f \-f-J 


( 26 ) 


(b) Point unit conductance on front portion of cylinder 
at any angle q> from 0° to 90°. The point of separation ' i s 
assumed to he located at <p = 90°. (See reference 24.) 



\ . il 3 ^ 

\ D / V 

, 90| / 


(27) 


Average Unit Conductances 
(a) Single Cylinder* 

f c = 0.211 T f 0 * 4 3 (^V)°- (28) 

D° « 4 


♦Based on reference 16, p. 222 of 2d ed. for a range of 
Reynolds numbers of 1000 to 50,000. 
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Figure 14.- Typical in-line tube bank. Number of rows, 4. 

Maximum G 0 based on minimum distance between 
tubes as shown in the figure. 



♦Based on the data in reference 8 for a Reynolds num- 
ber (based on the tube diameter) of 20,000. For Re of 
20,000 or greater, *he tube arrangement apparently has little 
effect on f c and equation (28) will yield results well 

within 10 percent of those in reference 8. For Re less 
than 20,000, however, the tube arrangement becomes more, im- 
portant, f c becoming lower as Sj, is decreased for a 

fixed value of Sf In most aircraft heater designs Re 
is in the vicinity of 20,000. If a design is contemplated 
with Re less than 15,000, p. 590 of reference 8 should be 
consulted. Note that the definition of F a in reference 8 

differs from that in this report. 
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TABL3 I.- TUBE ARRANGEMENT MODULUS FOR IN-LINE TUBE BANKS 


Number 
of rows 


1 


2 3 4 


5 6 7 8 


10 or 
more 


F q 1,00 1.10 1.17 1.24 1.29 1.34 1.37 1,40 1.42 1.43 

3 . 


(c) Banks of Tubes* (staggered.) 



1.25< j-<3 
5L 

1 . 25 < < 3 


Figure 15.-. Typical staggered tube baruc. Number of rows, 4. 

Maximum G 0 based on minimum distance between tubes. 
The position of the minimum distance varies with tube bank 
geometry. 


f„ = 14.5 X 10 
c 


-4 


A f 


0.43 & 0 


0 * 6 


, 0.4 


( 29a) 


“•‘Based, on the data in reference 8, for a Reynolds num- 
ber of 20,000. 
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TABLE II.- TUBE ARRANGEMENT MODULUS FOR STAGGERED TUBE BANKS 


Number 
of rows 


3 4 5 


6 7 8 


g 10 ° r 

more 


F n 1.00 1.11 1.23 1.31 1.39 1.45 1,48 1.51 1,53 1,54 

cL 


In' the foregoing equations. 


f 


c 0 


unit thermal conductance at stagnation point of cylin- 
der, Btu/hr ft 3 °F 



point unit conductance at any angle <P from stagnation 
point (up to point of separation, 0 < cp < 90°), 
Btu/hr ft 3 °F 


f c average unit conductance for any number of rows of cyl- 
inders, Btu/hr ft 3 °F 

3P angle from forward stagnation point, degrees. (*P must 
he less than 90°. ) 


If arithmetic average of tube wall absolute temperature 
and mixed-mean absolute air temperature, °R 

Ugc velocity in free air stream, ft/sec 

V density of air at temperature T^, lb/ft 3 

D outer diameter of tube, ft 

G 0 maximum weight rate per unit area for flow past tube 

banks (V/A), lb/ft 3 hr 

‘4 weight rate of air, lb/hr 

A minimum free area in flow path of air through tube 
banks, ft 2 


F a "tube arrangement modulus" for in-line and staggered 

tube banks , which is .mainly a function of the number 
of banks of tubes in the direction of air flow, 

(See table below each equation.) The magnitude of 
the modulus F a reveals the increase in the unit 
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thermal conductance for multi-row hanks with reference 
to single-row hanks, due to the increased turbulence in 
the • air stream. 


Exampl e 


What is the average unit conductance for the following 
hank of tubes? 



Minimum free area = (4 X 12) - ^4 x x 12^ = 24 sq in. 


= 0.168 ft 2 
3000 


0 o = 
u 0.168 


= 18,000 lh /hr ft 3 


J'a = 1.39 (table II, for staggered tube banks) 


Average "film" temperature = Tf = ~ + 460 = 591° 

D = 0.5 in. = 0.0417 ft 

Thus, the average unit conductance for the tube banks is! 


B 
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14.5 x 10 -4 x F a x T f ° * 43 x 
14.5 x 10' 4 x 1.39 x (591)° 
40.0 Btu/hr ft 2 °F 


D o • 4 

43 ( 1 8 OOP ) 0 * 6 

(0. 0417 ) 0 * 4 


' E . FORCED CONVECTION ALONG AIRFOIL SHAPES 

General 


As air flows over an airfoil, a laminar Boundary layer 
Is Initiated at the stagnation point and increases in thick- 
ness as the fluid proceeds. (See reference l4, p. 466.) At 
some point along the airfoil, depending largely upon the 
turbulence in the free air stream and on the airfoil design 
and surface finish, the laminar boundary layer will change 
Into a turbulent boundary layer. (See references l4, p. 482; 
35 and 36 . ) The growth of both the laminar and turbulent 
boundary layers is greatly affected by the pressure gradient 
along the airfoil. Except for the uncertainty concerning 
the point of transition between the laminar and turbulent 
boundary layers. It is a relatively simple task to predict 
the boundary-layer behavior around the airfoil. (See refer- 
ences l4, p. 156 ; 37* mad 38 .) An analytical prediction of 
the unit thermal conductance at a surface along which there 
exists a variable pressure gradient is, however, extremely 
difficult. See references 39 , 4o, and 4l for presentation 
of several analytical techniques.* An approximate solution 
is suggested in the following paragraphs. (See reference 

24. ) 


Along Leading Edge 

The b oundary- 1 ay er behavior near the stagnation point 
approaches that which occurs around a cylinder with a radius 
equal to the "equivalent radius** of curvature" of the air- 
foil leading edge. There is a difference in the behavior of 

*A comparison of these techniques is presented in a 
summary report soon to be published, 

**See example under sec, E, 



56 


/ 


NACA ARR No. 5A06 


the boundary layer over the actual airfoil as compared to 
that over the equivalent cylinder, of course, owing to the 
difference in the pressure distributions in fche two cases. 
As a first approximation, however, the leading edge of an 
airfoil may he considered as the front half of a right cir- 
cular, cylinder and the equations already presented for this 
case may he applied directly to evaluate the unit thermal 
conductance at the airfoil leading edge. (See ot. I. sec. 
D. ) 


Along Remainder of Airfoil 

As a first approximation the remainder of the airfoil, 
other than the leading edge, may he considered as an equiv- 
alent flat plate. The velocity along this equivalent flat 
plate is assumed to vary in the same manner as the velocity 
at the edge of the boundary layer of the airfoil. This ve- 
locity variation is easily obtained from the pressure dis- 
tribution about the airfoil .* 11 The origin of the equivalent 

*In order to obtain the velocity Uj at the edge of 

the boundary layer the following equations are utilized (ref 
erenee 42): 

From the Bernoulli equation, neglecting difference*, in 
elevation, 

*1 + 1 « *» + lissl • . (30) 

Y 2g Y 2g 

where 

P a static pressure at any point on the airfoil, lb/ft 
Y density of air in free air stream, lb / cu ft 
Uj velocity at edge of boundary layer, ft/sec 
P^, free stream static pressure, lb/ft 8 
u a free-stream velocity, ft/sec 
g gravitational force per unit 
Thus 

Hi = Ho, 


'1 - 


2g(Pi- - Poo) 


2 v 
u Y 

CO 


(31) 
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flat plate is assumed to be the stagnation point of the air- 
foil. The unit thermal conductance at any point x feet 
from the stagnation point then may be calculated from the 
equations for the point unit conductance along a flat plate 
presented in an earlier section. (See pt. I, sec. B. ) The 
exact point of transition from a laminar boundary layer to 
a turbulent boundary layer still remains in doubt. 


Recapitulation 


As shown in figure 17, in order to determine approxi- 
mately the variation of unit conductance over an actual air- 
foil, an equivalent airfoil made up of a cylinder and flat 
plate may be visualized. Application to the equivalent air- 
foil of the equations for the variation of unit conductance 
over cylinders and along flat plates then will. yield the 
desired estimate of the point unit conductance. The actual 
velocity variations over the airfoil, based on the pressure 
distribution, are to be utilized in calculating the point 
unit conductance! * 

Summary of Equations 

(a) Leading Edge (angle cp measured from stagnation point) 


f 


°CP 


0.194 T f 


o 


"(¥) 


0 




(33) 


(b) Laminap Boundary Layer (beyond leading edge) (x = 0 
at stagnation point) 


f Cx = 0.0563 T f °- 50 (^)°’ 


50 


(33) 


(c) Turbulent Boundary Layer (beyond leading edge) 
(x = 0 at stagnation point) 


f c = 0.51 T f 
x 


'.3 (u^) 


0 # 8 


O. 2 


(34) 


‘See preceding footnote. 
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where 

f c point unit conductance at leading edge of airfoil 
CP (0 < tp < 90°) Btu/hr ft 2 °F 

Tj arithmetic average of temperature of. airfoil surface 
and free stream air temperature, °R 

u^ free stream air velocity, ft/sec 

Y weight density of air at temperature Tf and pressure 
at given altitude, lb/ft 3 

D diameter of "equivalent" cylinder at leading edge of 
airfoil, ft 

f Cx point unit thermal conductance at point x "beyond 
leading edge, Btu/hr ft 2 °F 

u i velocity at edge of boundary layer at any point x 

along airfoil (obtained from static pressure dis- 
tribution over airfoil), ft/sec 

x distance along airfoil surface measured from stagnation 
point, ft 

Example 


Calculate the point unit conductance along an NACA 
S3012 airfoil with a chord of 6 feet for the following con- 
ditions? 


True airspeed 300 mph 

Altitude 15,000 ft 

Wing surface temperature .... 32° F 

Air temperature . 5° F 

Angle of attack 8° 


A sketch of the airfoil is shown in figure 17. The equiva- 
lent airfoil also is shown in this figure. By eye, the 
equivalent cylinder which best fitted the leading edge was 
found to have a diameter of approximately 0.282 foot. The 
f oil owing . equat i on now may be applied: 


(a) Leading Edge 

% - 0.194 (^iy- 5 [l 



(a ) 
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(b) Laminar Boundary Layer 

f Cx = 0. 0552 T f 0 * 5 0 


o. 5 


(c) Turbulent Boundary Layer 


f Cj =,0.51 Tf 0,1 <U »V. S 
5 + 32 


0.8 


The mean temperature = - — - + 460 = 478 R 


00 


(<?) 


The air density Y at 15,000-foot altitude is approximately 
0.048 pounds per cubic foot. (See appendix B.) The veloc- 


ity 


equals 300 miles per 

hour equal 

s 440 feet per 

second 

Thus 

the 

foregoing equations 

reduce to 




fc 

cp 

■ 34 - 4 i 1 - i&n 


(Leading Edge) 

(d) 


f C* 

/ \V 

= 0. 270 ( 2 

\ X / 

(Laminar 

Boundary Layer) 

(e ) 


f °x 


(Turbulent 

Boundary Layer) 

(f) 


In order to eval\iate the velocity u lt at the edge of the 
boundary layer, equation (30) is utilized. Thus: 


“l 



Poo) 


(g) 


The pressure distribution about the NACA 23012 airfoil, ob- 
tained from reference 43, is shown In figure 18. From this 

curve and equation (g) the velocity Uj at any point x. on 

the lower and upper surface of the airfoil can be readily 
determined. Substitution of these magnitudes of Uj in 
equations (e) and (f) allows the prediction of f c at any 

x for either the laminar or turbulent boundary layers, Tho 
resulting values of f c are plotted in figure 19. 





2.0 JO 40 SJO 6.0 7.0 

O/sZa/zce /?v/7? Sf&p/ioZ/oo /a//yf A/orrp 
A /rZa/V •J&r/ace - /*eZ: 


/v<?. /Sr 00/0 Sor A/A.CA. *230/2 A/rZb/V. 
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fig. /& - UNIT THERMAL CONDUCTANCE ALONG AN AIRFOIL 
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The point of transition between laminar and turbulent 
boundary layers is not well known, and this lack of knowledge 
limits the effectiveness of the prediction. In figure 19 
the point of transition was arbitrarily chosen at x =0.5 
foot on both the upper and lower surfaces in order to show 
how f<i may vary under certain conditions. In any actual 
design, information concerning the location of transition 
must be available in order to establish the variation of 
,f c with x. 


F. FINKED SURFACES 
General 


Fins, placed on a heat exchange surface, increase the 
rate of heat transfer by increasing the "effective” heat 
transfer area. When fins are used, they should be placed on 
the side of the exchange surface on which the thermal re- 
sistance between the fluid and, the surface is highest. Fins 
will have a negligible effect on the over-all conductance if 
placed on the side of the heat exchange surface having the 
lowest thermal resistance. (See reference 44.) 


Simplified equations for the performance of fins are 
shown. The general equation for the effective total conduct- 
ance, of the finned surface is (reference 45): 




n 



tanh 




(35) 


where 

(fA) e effective total conductance for finned surface, 

(Btu/hr)/°F temperature difference between base of 
fin and fluid flowing over firis 


n mimber of fins 

fji unit thermal conductance along fin, Btu/hr ft • F. 

, The magnitude of unit conductance fj? for fins 

may be • determined by use of 'data-in part I, sec- 
tions A, B, C, D, B, and G of this report. When 

radiant and convective heat transfer occur simul- 
taneously to the fin, fj = (f c + f r ) (Pt. I, 
sec, A . ) 
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L length of fin projecting into fluid stream, ft 

A cross-sectional area of fin perpendicular to direction of 
heat flow along fin, fts 

P perimeter of fin measured parallel to base of fin, ft 

k thermal conductivity of fin material, Btu/hr ft 2 (°P/ft) 

f unit thermal conductance along unfinned surface, 

Btu/hr ft 3 °F. The unit thermal conductance f u may 
"be determined by use of the data in part I, sections 
A, E, C, D, E, and G of this report. 

A u surface area not covered by fins, ft 2 

c 

Three particular forms are often met in practice; 



Fin perimeter, 

P = 2(s+l) ~ 21 

Fin cross-sectional 
area, A = si 


Figure 20.- Bectangular fins. 
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•Thus, from equation (35): 


+ 

where 


(fA), 


Tr'nD 


JVk f> 


4f F L 


tanh 


kD 


15 outer diameter of fins, ft 

3. Annular Fins (reference 46) 



(37) 


(fA) e = 


-nD c n v /2f F ks ^1 + j tanh 


ks 




(38) 


where 

outer diameter of cylinder to which fins are attached', ft 
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Example 

Determine the effective total conductance (fA) e for . 
tAe following aluminum finned surface. 


450 

lb/hr 










! 






t i! 

JJ 

5" 

/ rm rT7T 

rrrr 

/ / f t / / * / 

i nr 

n~rr 

n n 

TTTTTTT 7 

/ A"// // 

r f / f 7 '> f - 

r TTT 

rr* / / / i / ////'/// f 

r TT 


fyj * * f f f ' * ' ' r • t • '•*/// * f f j ft ’ r r / r r / 

Temperature of fin base = 800 °F 



(f.O e = 


nrrD 


•J Dkf p t anh 


/ 4 . f F L 

kD J 


+ Vu 


Inspection of the figure reveals that a reasonable procedure 
for the evaluation of fp is to consider the fins as a bank 

of tubes. The equations presented in part I, section D for 
tube banks then may be utilized to predict an average f'p. 


fp = 14.5 x 10" 4 F a T f ° 


4 3 


r 0.6 



r, o. 4 


' From table I 


F 


a ~ 


tv 


1.40 
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The film temperature, Tf = — — — — — + 460 = 1060° R (neg- 

o 

lecting temperature drop along fin as a first appr oximat i on * ). 
The minimum cross-sectional area for air flow is: 


A = 1,5 X 3 - 4 . g - P-V? . g . - X -i = 0.0243 ft 2 
144 

Maximum weight rate per unit area 


0. . -A 5 -°- = 18,500 


0 0.0243 


hr ft' 


D = 0, 25 in. = 0. 0208 ft 


Thus 

fp = 14.5 X 10“ 4 X 1.40 X (1060') °* 43 X = 69.4 Bt g a - 

( 0 . 0208 ) 0 * 4 hr ft F 


The exact evaluation of f u , the unit thermal conductance 

along the unfinned area, is difficult. As a first approxima- 
tion the f u may he calculated hy considering the unfinned 

surface as a section of a short tube (pt. I, sec.-C). 

Then : 


f 


a 


9.1 x 10" 4 T° 


3 



6N \ 
2 ) 


Fluid temperature = .T - 400 + 460 = 860° H 
Weight rate per unit area = G = 18,500 Ih/hr ft 2 
Tube length = l - 0.50 ft 


f = 9 1 x 10" 4 <fi6Q) 0,5 ( 1 S 5 00 ) ° * 8 
u * ‘ ( 0 , 5 0 ) 0 • 3 

The unfinned area, A^, = — - — — - — X 32 x 

u 144 4 


20.7 

0.35 2 

144 


Btu 

hr ft 3 °F 
= 0.114 ft 


Number of fins, n= 32 

*See reference 12, ch. II and reference 16, p. 232 of 2d 
ed. for temperature distribution along fin. 
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Diameter of fins D = 0,25 in. = 0.0208 ft 
Thermal conductivity of fin material, k = 133 BtU 

Length of fin = 1 in. = 0.0833 ft 


hr ft3 (D 


The term 


V kD 


' 4 x 69.4 x 0. 0833' 
133 X 0. 0208 


= 0.835 (dimensionless) 


"hen 


- j ~33«3x 3 . 


14 x 0.0208 


(fA)e~ 1 2 

= 14.5 tanh 0.835 + 2.36 


A 


0208 x 133 x 69 


'.4 tanh 0. 


835+ 20. 7x 0.114 


= 9.90 + 2.36 = 12,26 Btu/hr °F 


The fins, in this example, increase the effectiveness of the 
heat transfer surface appreciably. The heat transfer from 
the ends and radiant heat transfer between the fins are neg- 
lected here. A method of accounting for the heat transfer 
from the fin ends is presented in reference 47. 

The heat transfer rate, as computed herein, is slightly 
overestimated because the mean velocity is not uniform across 

ref erences 4g and 4 9 for experimental data 
♦ transfer, and velocity and temperature dis- 
tribution between longitudinal fins. 


of heat flow from the finned surface 


is de- 


If the rate 

sired, the value of (-fA) e must be multiplied ' by the mean 
temperature difference between the fin base- and the air flow- 
ing over the fins. The temperature distribution along the 
im need, not be considered. 


G. RADIATION 
General 


m an earlier section 
was pointed out that heat t 
companies convective heat t 
in parallel. The radiant h 
the fourth power of the abs 
body, becomes of great* impo 
bodies at high temperatures 
designs therefore, since th 
temperatures in the neighbo 
transfer may play an import 


of this report (pt. I, sec, A) it 
ransfer by radiation usually ac- 
ransfer , the two processes acting 
eat transfer, since it varies with 
olute temperature of the radiating 
rtance when heat transfer from 
is being studied. In some heater 
e heat transfer surfaces may attai 
rhood of 1000° F, radiant. heat 
ant role. 
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The heat transfer by radiation between two surfaces with 
temperatures Tj and T 2 may be calculated from equation 
(11) 

1r - 0-173 A r P AE [ (Jfc) - (Jfc) ] (35) 

The modulus ^aE» as discussed under part I, section A, 

modifies the equation for radiation -between Planckian radia- 
tors (black bodies) to account for the emissivities and ge- 
ometry of the radiating surfaces. The modulus a 

function of the areas of the surfaces ( A ^ and A 2 ) , the cor- 
responding emissivities e 1 and € 2 of the two surfaces 
and the relative geometry of the system. Values of Fas 

for the most common systems met in practice are given in 
table III. Other values for more complex systems may be 
found in references 16, pp, 54-60 of 1st ed. ; 1.9; and 20. 


TABLE III.- VALUES OF MODULUS ? AE 


System 


Modulus 


F A3 


Area to be used 
in equation (S) 


1. Two infinite parallel 
plates 


J^ + A 


- i 


A r = A, 


or As 


2. Completely enclosed 
small body (l refers 
to enclosed body) 


A r = Aj. 


3, Concentric spheres or 
infinite concentric 
cylinders (l refers to 
enclosed body) 


A ♦ (± - i) 

€ 1 ^2 V e O J 


A r = k x 


In table III, 

Aj surface area of smaller body, ft 2 


» 2 
A 2 surface area of surrounding body, ft 

4 

€ emissivity of smaller body (dimensionless) 

€ g emi.ssivity of surrounding body (dimensionless) 
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Values of surface emiaslvities € for various materials may 
be found In references 16 , pp. 5^-60 of 1st ed.; 19, 50, and 

51 . 

Inequation (8), Aj and A 3 are the surface areas of 
the smaller and larger bodies, respectively, as long as the 
surfaces do not contain re-entrant angles. In cases where 
the surfaces contain re-entrant angles, the area A, or A 2 
no longer signifies the surface area, and equation (8) is 
not exactly applicable. As an approximation, however, an 
effective "projected” area can be' used for A x or A 2 , For 
example, if the radiation between the surfaces in figure 24 
is required., the effective projected area can be used in 
equation (8) rather than the actual area of the inner surface 
The effective emissivity of the surface A r is higher than 

the actual emissivity of the metal because of the indentation 
For an exact treatment of this case see reference 19. 



Figure 24.- Effective area for radiation. 


In some designs, the space between the radiant surfaces con- 
tains substances such as water vapor, C0 3 , hydrocarbon 
vapors or other athermanous gases. These gases absorb and 
emit radiant energy in certain wavelength bands and thus the 
net exchange of heat between the surfaces is altered because 
of the characteristics of these constituents, An estimation 
of this effect may be obtained by use of the data presented 
in reference 11. 
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If either of the radiant surfaces is also transferring 
heat by convection, equation (12b) containing the equivalent 
unit conductance for radiation f r must be used. 


= (f c + f r ) A (t 1 - T a ) 


(12b) 


where 


0,173 A r r AB 


*F - 
* r 


- feT i 


A ( t ^ - T a ) 


(13) 


The calculation of the total heat transfer (convection 
plus radiation) is thus straight-forward when the surface 
temperatures T a and T 2 , and the gas temperature T are 

known. In most problems of a heater design or the prediction 
of its thermal performance, however, these temperatures are 
not all known. The magnitudes of T a and T 2 may be esti- 
mated by means of a heat balance on the particular surface - 
that is, by a consideration of the temperatures of the fluid 
on either side of the surface and the corresponding thermal 

resistances. Because these thermal resistances 1 

(to + fr) A 

are functions of the equivalent conductances for radiation 
f r which, in turn, are functions of the unknown temperatures 
T a and T 2 , the determination of these temperatures must 
be performed by trial and error. When these temperatures are 
found, the total heat rate may be calculated. 

In most heaters the radiant heat transfer will be found 
to be small compared to the convective transfer. In some 
instances, however, by special design, the heat transfer 
rate of an exhaust gas and air heat exchanger may be increased 
appreciably through the use of "irradiated convectors," which 
are placed in the ventilating air stream in "view" . of the 
surfaces heated by the hot gases, All the heat absorbed 
through radiation by theso surfaces is then transferred to 
the air by convection. The following example illustrates 
their effectiveness. 

Example 

For the system pictured in figure 25, where an irradiated 
convector is placed between the heated surfaces, which are 
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sides of two exhaust gas passages of an exhaust gas to air 
heat exchanger, what percent increase in heat transfer rate 
will result by the use of the irradiated convector? 


For this problem: let Tg = 1500° F and T a = 200° F, 

the absorption of the radiation between surfaces 1 and 2 by 

any intervening athermanous gases is neglected, the fluid 

rates and dimensions of passages are such as to yield the 

unit thermal conductances* f c = 20 and f c = 15 Btu/hr 

g a 

ft 2 °F For this case the modulus F A g = — ^ , The 

JL + -2 1 


e i 


e 


2 


emissivity of surface 1 will be taken to be 0.8 and that of 
surface 2 to be 0.9, so that F^g has the magnitude 0,735. 


-Irradiated convector 



*The values of f c may be obtained from the equations 
in pt , I , sec. B. 
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A heat rate balance on surface 1 gives the equation: 

Heat gained by surface 1 f r ora hot gas = Heat lost by surface 
1 to air by convection and to plate 2 by radiation. 


q T = f c A (T - t a ) = (f c + f r ) A (ti - T a ) ( a ) 

g ^ Si 1 


where 


f r - 



Yt, v , 

( T 2 Y ' 

A r F AE 

-MOO/ 

\100/ - 

A 

- T a> 


sur face 2 is? 



(b) 


Heat gained by surface 2 by radiation from surface 1 = Heat 
lost by surface 1 to air by convection. 


That is, 


f rg A (t 3 



f c A (t s 
c a 


T a). 



(c) 


Wh ere 



(d) 


T = t + 460° and, for this case, A r = A 

In equations (a), ( b ) , and (c) all the variables, except the 

surface temperatures T 2 and T 2 , are known. The simulta- 
neous solution of these expressions, therefore, will yield 
Tu and . However, the form of the equations is such as 
to necessitate a trial-and-error solution. As a first ap- 
proximation, a value of T x may be found from equation (a) 
by setting f r ~ 0, This magnitude would be the surface 
temperature if the radiant heat rate were 2 ero. This first 
value of T a may be substituted* into equations (c) and (d) 
in order to find an approximate value of Tg, With these 
first values of T 1 and T s the first approximation to the 
unit thermal conductance for radiation t* may be found 

from equation (b ) . When this value of f r ^ is substituted 

into equation (a), a better value of T a then can be calcu- 
lated. As before, a new value of T 2 may be found from 
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equations (c) and (d). By repeating this procedure, values 
of Tj and T 2 which satisfy the equations (a), (b), (c), 
and (d) are determined. The temperatures were found to he, 
tj = 855° F and t 2 = 410° F. Thus: 



0.173 A X 0.735 


/8 55 + 4 6 0\ 4 

/41 0 + 46 0\ 4 ' 

V 100 ) 

V 100 ) 


A (855 - 200) 


= 4.7 


Btu 
2 0 

hr ft 2 F 


The over-all unit thermal conductance when the irradiated 
convector is inserted between the two heated surfaces would he 

1 1 1 

(^ A )rad f Og A ^ c a + f r 

The corresponding value without the radiation plate 
present is 



The ratio of these two over-all conductances is 


(DA >rad . f ° e A f °. A 30 ^ fj 

UA _1_ 4. 1 + 1 

f °« A ( f c a + f r>) A 20 16 + 4,7 


The addition of the irradiated convector thus increased the 
heat rate by 17 percent. 


The reduction of the hydraulic diameter for fluid flow 
D in this case would have approximately doubled the static 
pressure drop due to skin friction. 
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II. SINGLE PASS HEAT EXCHANGERS 


The previous section presented the basic equation for 
the determination of the unit thermal conductances in heat 
exchangers utilizing as working fluids air and products of 
combustion of hydrocarbon fuels, The following section com 
bines these basic equations for the analysis and design of 
such heat exchangers. 


A. MEAN TEMPERATURE DIFFERENCE AND 
HEAT EXCHANGER EFFECTIVENESS 


Single-pass heat exchangers usually may be classified 
according to the mode of flow of the and cold fluids 

passing through them, as follows: 

1. Parallelflow . in which the two fluids flowing along 

the heat transfer surface move in the same direc- 
tion. 

2. Conntfirf 1 ow , in which the two fluids flowing along 

the heat transfer surface move in opposite direc- 
tion. 

3. Crossflow , in which the two fluids flowing along 

the heat transfer surface, move at right angles to 
each other. Two cases of this type of heat ex- 
changer are presented in this report. , In the 
first case each fluid is unmixed as it passes 
through the exchanger, and therefore the tempera- 
ture of the fluid leaving the heater section is 
not uniform, being, hotter on one side than on the 
other, (See fig. 28.) A flat-plate type heater 
(see fig. 36) approximates this type of exchanger. 
In the second case one of the fluids is unmixed 
and the other fluid is perfectly mixed as it 
flows through the exchanger. The temperature for 
the mixed fluid will be uniform across the section 
and will vary only in the direction of flow. An 
example of this type of heater is an unbaffled 
tube-bank crossflow type exchanger. (See fig. 

37.) In the example given, the air is mixed and 
the gas is unmixed. 
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Many exchangers utilized in practice do not fall exactly into 
any of the preceding classifications, "but are a combination 
of these; usually, however, one form of flow is sufficiently 
predominant to permit classification. 

The manner in which the temperature of the two fluids 
passing through the exchanger varies with distance along the 
heat transfer surface depends on whether the oxchanger is 
parall elf 1 ow , counterflow, or crossflow. Typical temperature 
variations for these flow types are shown in figures 26, 27, 
and 28 . 

It is evident that a constant temperature difference 
between the two gases does not exist. The effective tempera- 
ture difference which determines the rate of heat transfer 
between the two fluids is a function of the terminal tempera- 
ture differences (that is, the difference in temperatures be- 
tween the two fluids entering the exchanger (t„ - T„ ) and 

the difference in temperature between the two fluids leaving 
the exchanger ( T g 2 ” T a )), For the parallelfl ow and counter- 
flow exchangers the effective temperature difference is the 
familiar log-mean temperature difference. (See reference 12, 
ch. XIV; also figs. 26 and 27.) For the crossflow exchangers 
the effective temperature difference is a more complex func- 
tion of the terminal temperature differences.* (See refer- 
ences 1, p. 35 » 55 * and 56.) Charts for the determina- 

tion of the effective temperature differences for the three 
types of exchangers are shown in figures 29, 30, and 31. 

*The log-mean temperature differences and the effective 
temperature difference for crossflow are based on the follow- 
ing postulates : 

1, The over-all unit thermal conductance U must be 
uniform throughout the exchanger. However, as seen in pt. I, 
secs. H , C, D, and E, the unit conductances usually will vary 
somewhat bhroxighout the exchanger. Utilization of an average 
value of the over-all conductance usually allows the use of 
the log-mean or .effective temperature differences without ex- 
cessive error, (An exact evaluation of the rate of heat 
transfer is extremely complex when the conductances vary 
throughout the heat exchanger. Reference 52 presents an 
evaluation for the simple case in which the over-all con- 
ductance varies linearly with temperature. 

2. The relation between rate of heat transfer and tem- 
perature change of the fluid must be linear. This statement 
implies 

(Continued on p. 80 ) 
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PARALLEL PLOW 

(See fi?. 29- Also) 







T 0jSj 
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Cross Flow 



Fig. 26- Temperature Variation in Crossflow Exchanger 
Neither Fluid Mixed 
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Thermal Output of Exchanger 


The thermal output of a heat exchanger in. which no - 
is lost to the surroundings * may be expressed as: 

heat 



q. 

= (UA) AT mp 

Parallelf 1 ow 

(41) 

Form I 


q 

= (UA) AT mc 

Count erfl ow 

(42) 



q 

= (UA) AT mx 

Crossflow 

(43) 



q 

= *a c pa (T gl - T ai ) $ p 

Parallelf low 

(44) 

Form I I 


q 

= v a c p a ( T g 1 “ T aj ) ®c 

Count erfl ow 

(45) 

(Reference 

5?) 

q 

= W a c p a ^ T gi ’ T aj) $x 

Crossflow 

(46) 


.Form I will "be found convenient when all the terminal 
temperature differences are known, and is employed generally 
when designing an .exchanger to given specifications. Fig- 
ures 29,. 30, and 31 alloxir the graphical evaluation of the 
effective mean temperature differences. 

Form II is convenient when the temperatures of the 
fluids entering the exchangers and the over-all conductance 
UA are known, "but the temperatures of the fluids leaving 
the exchanger are not known. Form II is particularly 'useful 
in calculating the performance of a given heat exchanger at 
.various gas and air rates. -The two forms will yield exactly 
.the same results. Charts for the determination of the heatej 
effectiveness $> are shown in figures 33, 33, and 34. 

( Continued from p. 77 ) 

(a) An adiabatic exchanger or one which loses to the 

surroundings a constant, fraction of the heat 
transferred (reference 53) 

(b) Constant heat capacities ~ . 

(c) No change In phase (condensation, etc.) during any 

portion of the fluid's path through the ex-? 
changer. If the change in phase takes place 
over all of the fluid's path, however, the tem- 
perature usually will remain constant and the 
log-mean can' be utilized. 

•See references 53 and 57 for the case In which the 
heater loses heat to the' surroundings. 
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Tfi-U, 


1$ 30 . L - Mean Temperature Difference fur Counter Flow, 
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$P 



W<? <=p_ 

F/p. SZ ~ Meo-fer S/fccr/ve/?ess 


Air ft7/&//e/-/ r /otv fjrr/*7/7f*r. 
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In the equations for the heater output presented in 

equations (4l) to (46), the nomenclature is as follows! 

q thermal output of heater, 3tu/hr 

AT m p effective mean temperature difference for par- 

allelflow, shown in figure 29. For parallel- 
flow AT m p is the log -mean temperature dif- 
ference as defined in figure 26, °F 

AT ffiC effective mean temperature difference for 

counterflow, shown in figure 30, For 
counterflow AT mc is the log-mean tempera- 
ture .difference as defined in figure 27, °F 

AT mx , AtJ, x effective mean temperattire difference for cross- 

flow, shown in figure 31, For crossflow 
AT mx is a somewhat c ompl ex funct i on of 

terminal temperature differences, 0 F 
(See references 1, p: 35. an ^- 5^» 55» a &d 
5M 

T mixed-mean temperature of cold air entering ex- 

changer, °F 

T &2 mixed-mean temperature of hot air leaving ex- 

changer , °F 

Tg i mixed-mean temperature of hot gases entering 

exchanger, °F 

T - mixed-mean temperature of hot gases leaving ex- 

S2 changer, °F ’ ' . 

air rate, lb/hr 

Wg hot. gas rate, lh/hr 

Cp heat capacity of air at arithmetic mean mixed 

a temperature, Btu/lb °F 

Cp heat capacity of gas at ar i thme tic mean mixed 

temperature, Btu/lb °F 


effectiveness of paral 1 el f 1 ow heat exchanger 

T a 2 - Taj, 

that is, . = - temperature rise 

T £i ~ T a a 
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of cold fluid divided "by difference between hot 
gas and cold air temperatures at entrance to heat 
exchanger. (See fig. 32.) 

effectiveness of counterflow heat exchanger - that 
T — t 

is, $ = — — - temperature rise of cold 

C T - T 

Sl ai 

fluid divided by difference between hot gas and 
cold air temperatures at entrance to heat ex- 
changer. (See fig. 33.) 

effectiveness of crossflow heat exchanger - that is,, 

T - T * 

a 2 'aj 

d> x - ~ ~ temperature rise of cold fl,uid 

Si " 

divided by difference between hot gas and cold 
air temperatures at entrance to heat exchanger. 
(See fig. 34, ) 


over-all thermal conductance, (Btu/hr C F) which is 
defined by : 



UA ( f A ) a kA (f A ) g 


where 


(47) 


(fA) a , ( f A ) g total conductance on the air and 

hot gas side, respectively 
(Btu/hr °F) 

The total conductance equals the product of the 
unit conductance (f c + f r ) in Btu/hr ft 2 .°F, and 
the area of the heat transfer surface in ft 2 , 
for. unfinned heat exchangers.. For finned heat ex- 
changer, (fA) a , ( f A ) are the equivalent . t otal 

conductance on the air and gas sides, respectively 
(See example of finned exchanger calculation, pt. 
II, sec. 4 . ) 


thickness of heat transfer surface material measured 
in direction of heat flow, ft 

thermal conductivity of heat transfer surface 
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material, Btu/hr ft 2 (— ] . (In the usual heat ex- 

\ft/ 

changer the term L/kA is negligible compared with 
the other ' two thermal resistances.) 


From the equation for the thermal output of the heater 
it is noted that two variables affect this output more than 
the others: namely, the effective mean .temperature difference 
and the over-all thermal conductance. The mean temperature 
difference is usually fixed by design conditions, and may be 
readily determined from figures 29, 30, and 31. Then, to 
obtain a given heater output, a magnitude of UA equal to 
q/ATjn must be provided for by the design. 

Exampl e : 

A crossflow type heater in which neither fluid is mixed 
i6 to be designed to raise the temperature of 3000,pounds of 
air per hour from 10° to 400° F. The temperature of the hot 
gases available for heating the air is 1600° F, and the hot 
gas rate is 600C pounds per hour. What must be the over-all 
conductance of the heater, assuming no heat loss from the 
heat exchanger to the surroundings? 

If no heat is lost to the surroundings, the heat gained 
by the air is lost by the hot. gas. Thus 


*a °p a < T a s ' V> * “g °p g < T gl - T ga > 


The heat capacity of air at an average temperature of 205 ] 

= 0,241 Btu/lb °F. (See appendix, pt, IV, sec. A,). The 
heat capacity of the hot gas* at an approximate temperature 
of 1500° F is 0,277 Btu/lb °F. Thus the temperature of the 
hot gas leaving the exchanger is 

V c p 

a <T„. - T.J 


g8 


= T 


gl 


W g C Pg 


a a 


aj 


= 1600 - 3000 x (400 - 10)'= 1430 c 

6000 x 0.277 


F 


*The heat capacity of exhaust gases may be calculated 
from the data'given in appendix A for the pure components of 
the. mixture. The value 0,277 was taken to be that for pure 
air for this exampl.e. 
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From figure 31a, the mean temperature difference AT ax may he 

readily obtained. The parameters necessary for the use of 
the curves are 


K = 


T a 3 " T ai 400 - 10 


T 6i - T a a 


1600 - 10 


= 0.245 


and 


T gi " T g 2 1600 - 1430 


T gi " T ai 


1600 - 10 


= 0.106 


T - T 

At the intersection of <±> y = 0,245 and -1 — — = 0.106 

T - T 

• AT 

AT 

the value of — = 0.820 is obtained. Thus the effec- 

T g. “ T a, •' 

tive mean temperature difference between the hot gases and 
the air is: . 

AT = 0. 820 x, (1600 - 10) = 1300° F 


The output of the heater is to.be 

1 " W a (T »» - V.) 


= 3000 x 0.241 x 390 = 282,000 Btu/hr 


Thus the necessary over-all conductance UA is given by 


q _ 282000 

AT " 1300 

mx 


217 


B tu. 
hr °F 


A large number of heaters may be designed all of which 
will have a given value of UA, since any of the variables 
which control the over-all "conductance can be adjusted at 
will. The complete design of a heat exchanger for a given 
application, however, involves a series of compromises which 
can be made only by the designer in each special instance. 
Thus, in addition to the thermal output of the heater at 
given air and gas rates, items such as allowable maximum 
pressure drop through- the heater, manufacturing facilities 
and techniques, space, weight, life requirements, and so 
forth, have an important role in the final choice of a 
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heater design. A complete consideration of all these items 
cannot he presented in this part, hut methods for evaluating 
the thermal conductance UA for several given heater types - 
will he outlined. These examples will illustrate the appli- 
cation to heater design of the basic equations given in the 
first part of this report,* 

The thermal performance of most aircraft heat exchangers 
has been predicted within about 20 percent by means of the 
expressions contained in this report, and the predicted per- 
formances are -usually on the conservative side. Comparisons 
of measured and predicted heater performance for various 
types of heaters are presented in graphical form in refer- 
ences 29 , ^5, 58, 59, 60 , 61 , 62, 63, 64 , 65, and 66. Sea- 
sons for discrepancies between measured and predicted results 
are discussed in the texts of these reports. 


B. EXAMPLES 

Example 1 - "Fluted" Type Heat Exchanger 
Calculate the performance of the parallel flow "fluted" 


type 

heater shown i 

n figure 

35 

for the fo 

11 owing 

operat i 

ng 

cond 

it ions J 



• 





W a 

ventilating ai 

r rate = 

3000 lh /hr 





W 6 

hot gas rate = 

5000 lh 

/hr 






T gi 

temperature of 

hot gas 

69 

entering ex 

changer 

= 16 

00° 

F 

T ai 

temperature of 

cold ai 

r e 

ntering exc 

hanger 

= 10° 

F 



As a first app 

roximati 

on 

to the heat 

er perf 

orman 

c e , 

the 

actual geometry of 

the ends 

of 

the heater 

can be 

negl 

ec t 

ed 

and 

all the calculations ha 

s ed 

on the dim 

ens i ons 

at t 

he 

cen- 

ter 

of th,e heater a 

s shown 

in 

figure 35 f 

The le 

ngth 

of 

the 

heat 

er is measured. 

between 

the 

midpoints 

of the 

taper 

ed 

sec- 

t i on 

s , 









^Radiant heat 

t ransf er 

i s 

postulated 

to he 

zero 

in 

the 


four examples which follow. In actual heaters the radiant 
heat transfer is a small part of the total unless special 
irradiated convectors are used. (See example in pt. I, sec. 
G. ) 
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The flow systems for the air and gas sides then reduce 
to flow in a straight duct. 

The following pertinent dimensions can he evaluated from 
the data of figure 35 , 


L 

D H, 


length of duct = 1,1? ft 


4 A, 


hydraulic diameter on air side = — — 4 = 0.0577 ft 


% 


€ 


hydraulic diameter on gas side = 


4 A 


'€ 


* 0.0620 ft 


S 


A 


heat transfer area = 1,17 x 12.4 = 14.5 ft’ 


Rat i o 


L _ 1.17 

Du “ 0.0577 

a 


= 20.3 


The unit conductance for the heater on both the gas and 

air sides can be based on the data for flow in long ducts. 

(pt. I, sec. C),. (The Reynolds number for the air and gas 

flows will be found to be well over 10,000 and consequently 

equation (25) may be utilized to evaluate f c and f c , 

a g 

The equations are: 

For the air side 


* 5.4 x 10 


rp 0 • 0*8 

— 4 a 




0.8 


and for the gas side, 


f Cg = 5.4 x 10 


m 0 . 3p 0.8 

-4 


o.a 


H 


e 


|l + 1.1 (25) 

[l + 1.1 ] (25) 


The weight rates per unit area, G a , 0 g , may be readily cal- 
culated. The temperatures T' a and 


£ 


must be estimated 


and then checked with the final results , since the outlet 
gas and air temperatures are not known. 


• Reasonable values of T & and T g are as follows: 
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T a = 610 R 
Tg = 1990° R 

The weight rates per unit area are; 


G a = 7 V° 0 - °" = 15.300 — 

a 0-196 hr ft 2 

5000 iv 

G_ = = 23,700 — 12 

6 0.211 hr ft a 

Thus the air side unit conductance is: 


= 5.4 x 10- 4 X 610°° 3 X , 15300 !l! L + la x ^0577 1 

0. 0577 ° * 2 L 1.17 J 


= 15.3 Btu/hr ft 2 °F 


The gas side unit conductance is : 


f- = ,5.4 x 10 
c g 


-4 „ 1 99 0° * 3 x 237 00° * 8 


0. 0620 


0. 2 


\l + 1.1 X 0,0620 1 
I 1.17 J 

= 30.4 Btu/hr ft 3 °F 


The thermal resistance of the metallic wall is negligible. 
Thus: 


_ 1 _ 

UA 


A Vo. *c ■/ 14.5 \Z0.4 15.37 

a g 


or 


UA = 147 Btu/hr ft' 


In order to determine the performance of the heater, figure 
32 is utilized. The parameters necessary to determine the 
heater effectiveness 3> are: 1 

ir 



NACA ARR No. 5A06 


97 


B = 


UA 


147 


« W Q 0.241 x 3000 


= 0.203 


Pa a 


Wa - Cp a _ 3000 x 0.241 _ 0.533 

W - c-n 5000 X . 0. 277 

6 Pg 


From figure 32, the effectiveness <t> is 


$ p = 0.174 


Thus the heater output is: 

q = w a c p a ^ T gi “ T a a ) 

= 3000 x 0.241 (1600 - 10) 0,174 
= 200, 000 Itu/hr 


A .second, approximation can "be obtained, by using the 
computed temperature and recalculating the convective con- 
ductances. 

If the actual geometry of the ends of the heater has 
been considered in the calculation, due both to the change 
in heat transfer area and the variation in the unit conduct- 
ance- on the end surfaces, the heater output may be as much 
as 15 percent different than the amount calculated. To es- 
timate the temperature of the metallic surfaces, equation 
(16) (pt . I, sec. A) may be utilized. 

Jot test results on a fluted-type heater, see references 

29. 59. and 65. 


Example 2 - Flat-Plate Type Heater 

Calculate the performance of the crossflow "flat-plate" 
type heater in which both fluids are unmixed shown in figure 
36 for the following operating conditions: 
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W a ventilating air rate = 3000 lb/hr 
Wg‘ hot gas rate = 5000 lb/hr 

T gi temperature of hot gas entering exchanger = 1600° F 
,T &i temperature of air entering exchanger = 10° F 

As a first approximation to the heater performance,, the 
actual geometry of the ends of the heater can be neglected 
and all the calculations based on the dimensions at the cen- 
ter of the heater as shown in figure 36. The flow systems 
for the air and gas streams are similar to flow in straight 
ducts. 


The following dimensions may be evaluated from the data 
of figure 36. 

L a length of air duct = 0.583 ft 


Lg length of gas duct 'a 1.13 ft 


D H. 


4 A. 


hydraulic diameter of air duct = — = 0.0427 ft 


Du hydraulic diameter of gas duct 

g 


heat transfer area = 23,6 ft : 


a 
4 A 


g 


= 0.0516 ft 


S 


Ratio — — = 
D H 


0.583 
0. 0427 


= 13.7 (long tube, see pt. I, sec. C) 


The unit conductances for both the air and gas sides 
may be evaluated from the equation for flow in long ducts. 

The Reynolds number for both the air and the gas sides will 
be found to be between 5000 and 10,000, and therefore the 
equations in part I, section C are not exactly applicable 
but are used here because better equati ons "are not available. 
Thus, for the air side, from equation (25) 


= 5.4 x 10-4 


ip 0*3 q 0.8 
3 

d h ®* a 


F X 1 
l 1 + u r J 


0*2 


(25) 




(Mot fo aeo/e) 



100 
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o 

and for the gas side 



5 


4 x lCr 4 


m 0.3 p 0.8 

g g 



8 


X + 


1.1 


D 


_Hg 

L g 


(25) 


The weight rates per unit area can be readily calculated. 
The temperatures T a and Tg must be estimated and then 
checked with the final results, since the exit temperatures 
of the gas and the air are not known. Reasonable values of 
T Q and Tg are as- follows: 

T a = 610° R 

T - 1990° R 
& 

The weight rates per unit area are: 


G a . 


G g = 


300C 


19 x 0.0246 
5000 


= 6^20 


lb 


18 x 0.0158 


= 17,550 


hr ft 
lb 


•hr ft 3 


The air side unit conductance is: 


-4 w 610°* 3 X 6420° * 8 


f- = 5.4 x 10”' X 


0. 0427 0 * 3 


fi ♦ 1 . 1 x 0^37 1 
[ 0.583 J 


- 8.30 Btu/hr ft* °F 


The gas side unit conductance is: 


f , 5.4 x 10- 4 x 1990°- 3 x ns$0°- e + la x 0 

c g 0.0516°* 3 


fi + i.i x 9 . 5.11 ] 
l 1.13 J 

= 24.5 Btu/hr ft* °P 


The thermal resistance of the metallic wall is negligible. 
Thus 
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or 



— (— + -s--") 

23.6 \8 .30 24.5/ 


UA = 146 Btu/hr °F 


In order to determine the performance of the heater, figure 
34 is utilized. The parameters necessary to determine the 
heater effectiveness <I> X are: 


W 

V 


a 

S 


fPa 

c Pg 


3000 X 0.241 
5000 X 0.277 


0.522 

* 


UA 

w a c p a 


146 

3000 x 0. 241 


0,202 


From figure 34 the effectiveness $ x is; 


$ x =0.172 

Thus the heater output is ; 

q = W a c p a ( T gl ~ T ai ) = 3000 x 0,241 (1600 - 10) 0.172 

= 197,000 Btu/hr 


If the ends of the heater had been considered in the 
calculation, the heater output probably would be about 10 
percent higher than the value calculated. 

A great improvement in heater output can be obtained 
from the heater analyzed if the weight rate per unit area 
G a is increased. As noted in the calculations the thermal 
resistance in the air 'side is much greater than that on the 
gas side. Thus, a great improvement in heater output will 
result upon increasing the weight rate per unit area on the 
air side. This increase is readily accomplished by decreas- 
ing the cross section of the air passages. It is clear, ‘ 
however, that the decrease in air passage area increases the 
isothermal total pressure dro,p across the air side of the 



102 


NACA ARR No. 5A06 


heater. The proper size of the air passage is then a compro 
mise "between effective heat transfer and an allowable pres- 
sure drop. The metallic surface temperatures also should "be 
calculated in order to insure safe metal working tempera- 
tures (equation (16), pt. I, sec. A). 

See reference 66 for details of measurement and predic- 
tion of results in & flat-plate type heater. 


Example 3 - "Tube-Bank" Type Exchanger 


Calculate the thermal performance of the crossflow 
"tube-bank" type of heat exchanger shown in figure 37 for 
the following operating conditions. 

W a ventilating air rate = 3000 lb/hr 

Wg hot gas rate = 5000 lb/hr 

Tg temperature of hot gases entering exchanger = 1600° F 

T g a temperature of cold air entering exchanger = 10° F 

The flow system on the hot gas side consists of flow in 
straight ducts; on the air side the flow is over a bank of 
staggered tubes. 


The following pertinent dimensions 
from the data of figure 37. 

Lg length of tube through which the 

^Hg hydraulic diameter on gas side = 


Rat i o 




D 


H 


6 


13.8 


can be evaluated 

hot gas flows = 1 ft 
0,0783 ft 


D a outer diameter of tubes = 0.0833 ft 

A g heat transfer area on gas side = 9.82 ft 3 

t 

A a heat transfer area on air side = 10.5 ft 2 

Jj 

Since the ratio of — s. > 4.4 and the Reynolds number for the 
gas side = 19,900, equation (25) (pt. I, sec, C) may be used 
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Fig : 37 lube Bank Type Heat Exchanger 
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to calculate the unit thermal conductance on the gas side. 
The weight rate per unit area on the gas side is 

G = 5000 - 26,000 — !£-_ 

S 40 x 0.0481 hr ft 

The mean temperature of the hot gas is unknown, hut a reason 
able estimate is: 


T g = 1990° R 

(This value must he checked with the final calculated magni- 
tude . ) 


Then 


f„ = 5.4 X IQ" 4 -*£ 


T ® ® G 


'g 


D Hg ° 


' °» a / Dr \ 

A- 0 + l - 1 if ) 


f c = 5.4 X ID'* X 199°°^ a«000°" B ° ( x + x . x <k°ZH ) 
e 0 , 07 8 3 0 * 2 V 1 ' 

= 33.2 Btu/hr ft 2 °F 


The unit conductance for the air side may he calculated as 
follows* The Reynolds number for the flow over the tubes in 
the tube bank is found to be 15,800. Equation (29) presented 
in part I, section D is then applicable. 

r o.e 

4 0,43 ^Oo 

f<j a = 14.5 x 10- 4 F a Tf ft % 7 C ~ 

£t 

The tube bank has 10 rows of staggered tubes. Thus from 
table II, F a « 1.54. The minimum area of flow for the air 
is shown in figure 37 and equals? 

A = ( 1 * 7 * 12 = 0.325 ft 2 

\ 144 / 

Thus the maximum weight rate per unit area is 

G 0 = - 5 = 9230 lb/hr ft 8 

° a 0.325 
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The approximate surface temperature of the tubes can he pos 
tulated to he midway between the mean gas and air tempera- 
tures. Thus the film temperature Tf can he calculated as 
follows : 


T 


f = 


( 15 0 j_ 1 530 , 160 


i- + 460 

2 


955° R 


Then 


. . 14.5 > 10-“ » 1.54 X 9560,43 11 37.7 Btu/hr ft* 

' a 0.0833°* 


(Since f _ = f„ , the assumption that the tube surface 

- c g c a 

t emperature is the arithmetic mean of the air and the hot 
gas temperatures is justified.) (See equation (16), pt. I, 
sec. A. ) 

The thermal resistance of the metallic wall is negligi- 
ble, so that: 

1 _ 1 + 1 _ 1 + 1 

UA (f„A) ( f e A ) 27.7 x 10.5 32.2 x 9.82 

c a G g 

or 

UA = 153 Btu/hr °P 

The hot gas, which flows through the tubes, is unmixed; 
while the air passing over the tubes is mixed. The heater 
is therefore classified as a crossflow heater, gas side un- 
mixed, air side mixed, (See pt, II, sec. A, item 3.) In 
order to determine the performance of the heater, figure 34 
is utilised. The parameters necessary to determine the 
heater effectiveness are: 

B UA ,- = i^3___ , 0 Q 

c,, W. 0. 241 X 3000 
a 


°Pa _ 3000 x 0,241 _ 0 .522 

W_ c n ~ 5000 x 0.277 
S Pg 

i 

Prom figure 34 the effectiveness $ x is 
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= 0.180 

Thus the heater output is: 


= 3000 x 0.241 x (1600 - 10) x 0.180 = 207,000 Btu/hr 


on A £ f ” ea8 ? PeBents and P redi ction« Of results 

on a tubular-type heater is given in reference 64 . 


Example 4 - Finned-Type Exchanger 


Calculate the performance of the 
exchanger shown in figure 38 for the 
diti ons : 


crossflow finned-type 
following operating con- 


W a ventilating air rate = 3000 lb/hr 
w g gas rate = 5000 lh/hr 


T g 1 temperature of hot gases entering the exchanger = 1000° F 
T a a temperature of air entering the exchanger = 10° F 


The exchanger consists of an aluminum casting with longitu- 

wJth X ^f inS l n th ® center along which the hot gases flow, and 

flow* Zl en ^ finS ° n the 0utside * along which the air 

flows at right angles to the directi on of the hot gases. The 

hot gases flow through a system consisting of a straight 
duct. As a first approximation the air flowing along the 
circumferential fins can he considered as flow in a curved 
duct, and as a further approximation the equation for the 
unit thermal conductance in straight ducts utilized to cal- 
culate f c . The following pertinent dimensions now can be 


a 


g 


cross-sectional area for air flow (sec. A-A) = 

(12 X 1.45 - 40 X 0.14 X 1.25) (2) 

144 — 1 *■ 0.145 ft 3 

cross-sectional area for gas flow (sec, B-B) = 


it x 36 
4 X 144 


30 X 1.30 3 

144 * 16 * 


0.145, ft 3 



V<7 fms on Air S&e 
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P a wetted perimeter on air side (sec, A-A)- = 

(1,25 X 3 x 40 + 24 + 2.90) 2 _ 21 2 ft 
12 

Pg wetted perimeter on gas side (sec. B-B) = 

1,30 x 30 X 2 + rr x 6 - 30 x I + 30 x I 
4 8 

1 2 

Hydraulic diameter on air side 


= 7.76 ft 


„ 4 A- 4 X 0.145 

D «a ' -yf - 21. Z 

Hydraulic diameter on gas side 

4 A» 4 x 0.145 


= 0.0273 ft 


7.76 


= 0.0747 ft 


Length of duct on air side 


L * ■ ff fa - °- 916 ft 


Length of duct on gas side 


L_ = 1 ft 


Determination of unit" thermal conductances: 

The Reynolds number on the air side = 10,000 

The Reynolds number on the gas side = 23,700 

Thus equation (25) (pt. I, sec. C) is applicable for 
the determination of the unit conductance on both the gas 
and the air side: 


Air Side 


c a 


= 5.4 x 


4 T 

10“ ZSL 


0 , 3 


0 . 8 


DH ; 


O . 3 


( 


1 + 1.1 


J H. 


) 
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The mean temperature T a is not known, "but a reasona- 
ble estimate is 550° R. This estimated magnitude must be 
checked with the final calculated results. The rate of flow 
per unit cross-sectional area G a is! 


G a = -?-QP9. = 20,800 — 

a 0.145 hr ft 2 


Then 
f , 


*e 


5.4 X IQ* 4 x ^0°- 3 x 208 


0. 0273° 


will (i + i.i 

* 2 \ 0,92 J 


= 22.0 Btu/hr f t 


a o. 


This unit conductance can be assumed to exist along both the 
fins and the unflnned area. 

Gas Side 

From equation (25) 


*e 


= 5.4 X lO - * 


rp 0*3 p 0*8 


DH 


o. a 


1 + 1.1 


■'H, 


S 


The mean temperature Tg is not known, but a reasonable es- 
timate is 1420° R. This estimated magnitude must be checked 
with the final calculated results. The rate of flow per 
unit cross-sectional area G a is 


Then 

fc S 


G 


a 


5000 

0.145 


34,500 


lb 

hr ft 2 


5.4 x 10” 4 x 


1420° * 3 x 34500°* 6 
0. 07 47 0 * 3 


^1 + 1.1 


0. 0747^ 

1.0 J 


= 37.0 Btu/hr ft 2 °F 

This magnitude of the unit conductance can be assumed to ex- 
ist along both the fins and the unfinned area. 



NACA ARR No. 5A06 


I 10 


Effective Conductance of the Finned Surfaces 
Air Side 

The extended surface on the air side of the exchanger 
consists of circumferential fins. Thus equation (38) Opt. I, 
sec. F ) is applicable for the evaluation of the effective 
thermal conductance on the air side of the heat exchanger 
surface. 


(fA) ea = tt D c n^y^SfjMcs" ^1 + tanh 



t fn Ai 


u-°-u 


From figure 38 and the previous calculations: 

D c diameter of cylinder to which fins are attached = 0.531 ft 
n dumber of fins = 40 

fj> unit thermal conductance from fin to air = f c& 

- 22.0 Btu/hr ft 2 °F 

k thermal conductivity of fin material (aluminum) at an 

average temperature of 200 ° F - 120 Btu/hr ft 2 (°F/ft) 

s thickness of fin = 0.0117 ft 

L length of fin perpendicular to cylinder = 0,104 ft 

f a unit thermal conductance from surface of cylinder to air, 
approximately equal to f c = 22.0 Btu/hr ft 2 °E 

£L 

A u area of cylinder not covered by fins = rr X 0.531 x 1 

- 40 x 0.0117 = 1.20 ft 2 

Substituting the foregoing values in. the equation for (fA) ea 
yields: 

(fA)ea = tt x O .531 x 40 J 2 x 22.0 x 120 x 0.0117 ( 1 + ) tanh 0.582 

\ O. 531 / 


= 1.20 x 22.0 


( f A ) ea = 329 + 26 = 355 Btu/hr ft 2 °F 
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Gas Side 

On the gas side of the exchanger the extended surface 
consists of rectangular fins parallel to the direction of gas 
flow. Thus equation (36)- is applicable for the evaluation of 
(fA) e g. 


(fA) e g * nl J 2skfp tanh 


2f -a L‘ 


ks 


^u-^u 


From figure 38 . and previous calculations; 

L 0.108 ft 
n number of fins, 30 

l length of fins in. direction of gas flow„= 1 ft 

s thickness of fins = 0.0156 in. 

k thermal conductivity _ of^ aluminum at 600° F 
= 140 Btu/hr ft 8 

fp unit thermal conductance from gas. to fin * f c 


(If) 


= 37.0 Btu/hr ft 2 °F 


g 


f u unit thermal conductance from gas to part of surface not 


covered by fins = f =37,0 Btu/hr ft 


3 On 


'g 


A^ = area of cylinder not covered by, fins = n x 0.5 

. HIM.. . 0.96 ft 2 

13 ... 

Substituting the preceding magnitudes of the variables into 
the equation for (fA) e yields': 

- ‘ - ^ ’ r v;. * 3 ; ....... . , . . . ' ‘ . 


(fA) eg c= 30 x: 1 y 2 x 0 v 0l:5 6 x. 140 x 3 7. tan'h^O. 63 0 + Q. 95 x >37 : - 


(fA) e = 213 + 35 = : *248; Btu/hr '-^•F . : 


eg 


The over-all conductance of the heat exchanger may now be 
calculated: 
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-1— + -L- = 0. 00686 

355 248 


Thus 


UA = 146 Btu/hr F 


Because "both gases flow between fins and are thereby 
prevented from mixing, the heater is a crossflow type heater 
with neither fluid mixed. (See pt. II, sec. A, item 3.) In 
order to determine the performance of the heater figure 34 
is utilized. The parameters necessary to determine the 
heater effectiveness $ x are: 


Wc 

a Pa _ 3000 x 0.241 _ 0 g50 
W g c pg " 5000 x 0.263 ' 


_UA_ 
W a °Pa 


— = 0.202 

3000 X 0.241 y 


From figure 34 the effectiveness is: 

3> x = 0. 172 

Thus the heater output is: 


= W, 


'a 


(T 


€i 


~ T a a > 


q = 3000 X 0.241 (1000 - 10) X 0.172 
q = 123,000 Btu/hr 


It should he noted 'that although the ef f ectiveness of 

this heater Is as high as those determined In the other ex- 
amples, the heater output is smaller because of the lower 
temperature of the entering hot gas. For test data on finned 
exchangers, see references 58, 60,„ and 63 * 
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III.. HEATER PERFORMANCE IN FLIGHT 

A. HEAT REQUIREMENTS 
1. Cabin Heating 

The heat requirements for cabin heating depend upon the 
heat lose through the cabin walls and the rate of air leakage 
into the cabin. (See reference 67.) In order to establish 
-the heater requirements, a heat balance can be performed on 
the cabin of an airplane as shown in figure 39. 


Wg (hot gases) 



Figure 39,- Mass and thermal balance on airplane cabin. 

Then 5 


V g c Pg ( T gi “ T g 3 > + c P a ^ V a + W L) ( T a “ T c) ~ <lw = 0 

but 

W g c P g (T Si - T g 3 ) = qn (Heater output) 


Thus : 
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qjj - (W a + W L) c p a ( T c “ T a^ + 4w 


or 


4H 


T _ T 
c 1 a 




L c Pa + 


l w 


( T c - T a ) J 


+ V„ c 


a '-' 0 . 


since q w , the rate of heat flow through cabin walls equals 
q w = (UA) (T - T & ) ' the following relation may he written: 




( T - T ) 
c a ' 


= (necessary heater output in (Btu/hr) per 

degree difference in temperature between 
cabin air and outside air) 


’ W 1 c P a + UA ] + "'a ° Pa 


or 


= UA 1 + 




(W L + w ) 


UA 


'Pa 


“ 


(48) 

(4.9) 


where 


rate of flow of leakage air into cabin, lb/hr 

Cp unit heat capacity of air at constant pressure, Btu/ 
a lb °f 


UA over-tall thermal conductance between cabin air and out- 
side air, Btu/hr F 

W a ventilating air rate for a ram operated heater, lb/hr 

The multiplier of UA in equation (49) i 
heater output required when there is air 
through the heater W a or through cracks 

would be required if there -was no leakage 
magnitude of this multiplier is a measure 
of an airplane cabin. 

The magnitudes of UA can he readily estimated by uti- 
lising the equations presented in reference 67, or by install 
latlon of heat meters (references 68 and 69) on the cabin 
walls. Obviously UA can be greatly decreased by the use 
of insulation. 


s the ratio of the 
leakage (either 

Wj,) to that which 

of any kind, The 
of the "tightness" 
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The leakage air rate Wj, usually is unknown, hut fdr a 

given airplane may he determined hy measurement, for equation 
(48), By measuring in fligh-t all other terms in the equation 
except W^, this magnitude may he easily determined as a 
function of altitude, airplane speed, and - so forth. Efforts 
should he made to reduce the leakage term as much as possible 
in order to reduce the thermal output of the heater required 
to maintain a given cahin air temperature,* 

It is instructive to note that, particularly if the 
leakage air rate Wj, is made negligible, it is advantageous 
to design a heater with a low ventilating air rate W a , 
since the smaller W a , the, smaller the heater output re- 
quired^ If the leakage is large, however, decreasing W a 
has a small effect on the necessary heater output, 

2. Wing Anti- Icing 

In the case of wing anti-icing, the conductances over 
the airfoil and in the air ducts may he estimated hy means 
of the equations presented in this report. Valuable flight 
data on wing anti-icing systems will he found in the reports 
from the Ames laboratory. (See references 70, 71, 72 , 73 , 

74, and 75 *) A valuable report on wing anti-icing has been 
written recently hy Myron Trihus. (See reference 82.) 

B, CORRECTION OF THE PERFORMANCE OF HEAT EXCHANGERS 
TO ALTITUDE CONDITIONS 
1, Thermal Performance . 

The thermal output of a heater, as discussed in part II, 
may be written as: 

*If the temperatures in the cahin are sufficiently uni- 
form, it may he possible to de t ermine ' - b oth UA and W* L by 
direct flight measurement as follows: At one altitude and 

airplane speed measure W a , T c ,T a a t two different values 
of Q.H • This will allow equation (48) to he written twice 
with two unknowns, UA and W^, Since these variables are 
practically unchanged in the two tests, they may he solved 
for. By performing these tests at various altitudes and 
airplane speeds the variation of UA and W^ with these 
variables may he determined. 
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= W a c p 


(t 


a 


Si 


- T ) $ 

a l ' 


(50) 


where $ , the effectiveness of the heater, a function of 


( 


UA 




W a °P a / 


and 


V/ c ' 

a Pa 

w s c Pg, 


is shown in figures 32, 33, and 34, 


for parallelflow , counterflow, and crossflow exchangers. 
Inspection of equation (50) reveals that for fixed values 
of V a Cp^ and c-p^, the heater output is proportional 


W 6 C P* 


'a * ^g' 

to the ratio of the entrance temperature difference and . the 
heater effectiveness Thus, for fixed weight rates* 


^alt _ ^ T gi ~ T a^alt $alt 
*lab ^ T gl " T a a ^lab $lab 


The ratio of inlet temperature differences is self-explanatory. 
This ratio depends only on the temperatures existing during 
the laboratory test and during the actual operating condi- 
tions at altitude, The magnitude of the atmospheric pres- 
sure does not affect this ratio. 


The second ratio is more complex, hut, as shown in ref- 
erence 58, for equal weight rates of gas and air It is in- 
dependent of atmospheric pressure and slightly dependent 
upon the mean absolute temperature of the two gases. As 
shown in reference 58* for the maximum temperature variations 


met in practice, the ratio 


$alt 

does not' vary from unity 


more, than 10 percent, Because of the complexity of evalua- 
tion of this ratio, it is suggested that for a first esti- 
mate the heater output can he corrected to -any altitude with 
sufficient accuracy by multiplying by the temperature differ- 
ence ratio only. Thus for fixed gas and air. weight rates 


q al t _ ^ T g : ~ T a 1 ^alt 

q lab (T g x * T ai ^lab 


*Cp and Cp vary slightly with temperature and are 

practically independent of pressure in the range under con- 
sideration. 
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where 

<lalt 

<U ab 


heater output at any altitude for given values 
of air and gas rates, Btu/hr 

heater output obtained in laboratory (or calcu- 
lated), Btu/hr, for the same. values of W g , 
W a as required for q a ]_t 


(t -t ) difference between the temperature of hot gases 
1 1 a and cold air entering exchanger at any alti- 

tude, °F 


( T ei** Tft i >lat 


difference' in temperature between hot gases and 

cold air entering heater, in laboratory test, 

op 


If a more precise correction is required, reference 58 should 
be' consulted. 


2. Pressure Drop across Heater 


a. I s othermal , - The basic pressure drop measurement re- 
quired to establish heater performance is the isothermal 
total^pres sure drop across the heater. The isothermal total - 
pr e s sure drop represents the loss due to frictional forces 
such as skin friction, sudden expansion, sudden contraction, 
and go forth, and is an irrecoverable loss. Pressure drops 
due to the acceleration of the fluid may be recovered. 

The isothermal total-pressure drop may be obtained ex- 
perimentally in two ways ; 


fl) By traversing the duct before and after the heater 
with a total-head tube during isothermal flow'. 
The total-head loss (reference 7# p. 206 and 
reference 76) for a heater with circular inlet 
and outlet ducts then will be: 


dF„ -k = (p - PO- 

a-b a b 


3600 y 
2g 


/ 2tt v 3 rdr - / 2rr v^ 3 rdr 


W 


( 53 ) 
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The integrals indicated must be evaluated graphically, uti- 
lizing the pitot tube data to establish v a , v^ as a func- 
tion of r. An approximation of these integrals may be ob- 
tained from measurements with a pitot tube at appropriate 
points v across the pipe section. (See “ten-point method" in 
reference 77.) 

In- oquati on (53) 


^a-b 

Pa 

p b 

V 

g 


isothermal frictional pressure loss between sections 
a and b, lb/ft 2 

static pressure at section a, lb/ft 3 
static pressure at section b, lb/ft S 
density of fluid, lb/ft 3 


gravitational force per unit mass » 32.2 lb 


//lb sec 2 \ 

A—x~) 


velocity of fluid at any radius of pipe r at sec- 
tion a, ft/sec 


v-j) velocity of fluid, at any radius of pipe r at sec- 

tion b , f t /sec 


r any radius, ft 

r a inside radius of pipe at section a, ft 

r-b inside radius of pipe at section b, ft 

W weight rate of fluid, lb/hr 

(2) If the areas at a and b are equal and if the 
velocity distribution acr oss the two sections 
is postulated to be the same, equation (53) 
reduces to; 

AFa-b (friction pressure loss) 

= P a - Pb (static pressure drop) 


Thus , as an approximation, the isothermal friction pressure 
loss across a heater may be obtained by the static pressure 
drop measured at sections of equal areas. 
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to. ffonlsothermal . 


pressure 


- Once the isothermal frictional 
has been obtained for a series of 


weight rates, the nonisothermal static pressure drop for 
the same weight rates at any altitude may be readily calcu- 
lated toy means of the following equation (reference 16 , p. 
130 of 2d ed., and reference 78): 




n on-i s 0 




where 

P static pressure at section a, entrance to heater, 

a lb/ft 2 

Pv static pressure at section b, exit from heater, 

lb/f t a 

hP a _-^ isothermal friction pressure loss through heater at 

weight rate W^go* temperature Ti s 0 * and pressure 
Pi s0 , lb/ft 8 (This friction pressure loss in- 
cludes any irrecoverable losses from sudden con- 
traction, expansion, etc., as well as skin fric- 
tion, ) 

T ft temperature of fluid entering exchanger, °R 

temperature of fluid leaving exchanger, °R 

Ti g0 temperature of fluid during isothermal pressure drop 
determination, °R 

Piso average pressure in heater during isothermal pressure 

. drop determination, lb/ft 2 abs. 

P average pressure in heater at any altitude, lb/ft 2 abs. 


W^ g0 weight rate of fluid for which isothermal and non- 
isothermal pressure drops are being determined, 
lb/hr 
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R gas constant for air = 53.3, ft-lb/lb °R 

g gravitational force per unit mass = 32.2 lb 

constant cross-sectional area of portion of heater in 
which heat transfer takes place, ft 2 

A a cross-sectional area of inlet duct to heater at which 
static pressure P a is measured, ft® 

At) cross-sectional area of heater outlet duct at which 
static pressure is measured, ft 8 

The first term on the right of the e^ual sign 


lb sec 8 N \ 

ft J 


A P 


a-b 



represents the irrecoverable, non! so t hemal pressure loss 
due to friction, including sudden expansion, contraction, 
and so forth. This term, which also can be written 



is used to predict the non! sothermal frictional pressure 
losses from the isothermal values. 


The ratio 



is an approximate correction for 


the change in friction factor with change of Reynolds 

number (Reynolds number changes are caused by changes in 
absolute viscosity as the fluid is heated or cooled)* 


The ratio ffir) is an approximate correction for the 

variation of pressure drop with changes of density caused by 
temperature and altitude effects* Because of the fact that 
an arithmetic average density is used, the correction (for a 
fluid being heated) is too high for the heater-entrance con- 
traction losses and too low for the heater-exit expansion 
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Section 


Area 


Temperature 
(iso. ) 


Temperature 
(noniso. ) 


Kean velocity 


Static pressure 


Average absolute 
static pressure 


Specific volume 

(noniso . ) 



Figure 39a.- Isothermal and noniso thermal flow 
conditions in heater. 
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losses. These errors' 1 probably partially compensate one an- 
other. A slight error is made, also, when AF a _ b is multi- 
/ T V® *13 

plied by ( ) , because expansion and contraction 

\ T iso J 

losses should not be corrected for changes in viscosity with 
temperature , 

The seoond term to the right of the equal sign of equa- 
tion (54) represents a pressure, drop due to the acceleration 
of the fluid, which resul-ts both from changes in cross- 
sectional areas and changes in density due to heating or 
cooling. This term 


V 3 600 ) 


R T. 


\ 3 6 0 0 / 2g A n 2 P 
may be written as 




v 


(2 P * • I P * U a 3 ) + (s tt V - I P* u » a 3 ) 


or more simply 


( 1b - + (c *hv - <lh > 


where q is used to denote the velocity pressure p u /2 
Equation (54) may thus be rewritten as: 


f ^i s 0 \ / ^a v \ ° * 




ISO 


13 


+ ( Rh b “‘lh a ) (54a) 


♦These errors are reduced hy use of the following ex- 
pression: 



in which the pressure loss through the heater shown in fig. 
39a is divided into a contraction term, a friction term, and 
an expansion term. 
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The term (F a + q a ) - (P^ + q-j) ) is the difference in total 
pressures at points a and b as measured by an impact 
(pitot) tube. Since total pressure represents the energy 
available for pumping this fluid, equation (54a) probably 
has more physical significance than its equivalent, equation 
(54). It should be noted that the noni soihermal frictional 
(or irrecoverable) pressure loss between sections a and b 
(the term in equation (54a) involving &F a _t,) differs from 
the difference of total pressures by the term ( - q^ ) 

which is the change in velocity pressure through the heater 
resulting from the heating or cooling of the fluid. The 
term (q^ ~ *lh a ^ for a ^ r ^ e i n S heated represents a loss, 

which cannot be "recovered" except by cooling the fluid as 
it passes through the discharge duct.* It is not possible 
to recover the loss qj^ - q_h a by diffusers and other me- 


chanical means. Thus, if the discharge duct from .a heater 
is adiabatic, the term - q^ in effect represents an 

irrecoverable loss which should be charged against the 
heater. If the fluid is cooled in the discharge duct, how- 
ever, a gain in velocity pressure will result. which can be 
less than, equal to, or greater than the loss q h ^ - q h , 

depending on the amount, of cooling; that is; cooling the 
fluid is equivalent to introducing a pump in the discharge 
duct. In practice, however, the term q h - q. usually is 

n b n a 

small and can be neglected, It should t>e emphasized, however 
that equation (54a) reveals that in nonisothermal flow the 
frictional pressure loss is not exactly equal to the total- 
pressure difference. 




C. ALTITUDE PERFORMANCE OF HEATER AND DUCT SYSTEM 


The analysis of the performance of a ram (or fan) oper 
ated heater and duct system (see fig. 4o) as a function of 
true airplane speed and altitude is presented in detail in 
reference JS. The two basic equations necessary for the 
analysis are presented: 


*This is accomplished in a wing de-icing system. 



Figure 40.- Elements of heater and duct system. 
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where 


APf an total pressure rise across a fan which may he placed 
in inlet duct , lb/ft 2 

n exponent to account for friction loss variation with 

W due to expansion and contraction and to skin 
friction (its value will be between 1.75 and 2,00; 
nearer to 1.75 if the major loss is due to skin 
friction.* (May be obtained from a plot of Wj gQ 

versus isothermal frictional pressure loss through 
system. ) . ) 

Pi absolute static pressure in free air stream before 

air scoop, lb/ft 2 

*See references 15 . 29, 45, 58, 59. 60 , 6l, 6 ? and 63 
for typical values of the exponent n. 
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true air speed of airplane + air velocity pro- 
duced by propeller (velocity of air stream 
relative to airplane, ahead of air scoop), 
ft /sec 

specific volume of air in free air stream be- 
fore air scoop, ft 3 /lb 

gravitational force per unit mass 
= 32.2 lb/( x * f 3 t aC ~ ^ 

absolute static pressure at point of air dis- 
charge, lb/ft 2 

velocity of air relative to airplane at point 
of air discharge, ft/sec 

specific volume of air at point of air dis- 
charge, cu ft/lb 

total pressure in free air stream before scoop, 
lb/ft 

air flow through duct, lb/hr 

air flow through duct for which isothermal 
total head losses, ^ F i-r 2 . AF 2 _ 3 » AF 3 _ 4 , 
AF 4 „ 5i etc., were determined, lb/hr 

gas constant for air = 53.3 ft lb/lb °R 

absolute temperature of air in free air stream, 
°R 

absolute temperature of air, just inside air 
scoop , °R . 

absolute temperature of air at entrance to heat 
exchanger, °R 

cross-sectional area at section 3-3, entrance to 
heat exchanger,' ft 2 

mixed-mean absolute temperature of air at sec- 
tion 4-4, exit of heat exchanger, °R 
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t 4 mixed-mean temperature of air at section 4-4, exit of 

heat exchanger, °F 

A 4 cross-sectional area at section 4-4, exit of heat ex- 

changer, ft 3 

T 5 mixed-mean absolute temperature of air after passing 

through nonisothermal duct just before final dis- 
charge section, °R 

A g cross-sectional area of duct at section 5-5, just be- 

fore final discharge section, ft 2 

average static pressure in duct system during the 
isothermal total-pressure drop test, lb/ft 3 

frictional pressure loss between the free air stream 
and entrance to air scoop for isothermal condi- 
tions specified by P is0 , T is0 , W ia0i lb/ft 2 

frictional pressure loss between entrance of air 
scoop and entrance to heat exchanger, for iso- 
thermal^ condi ti ons specified by P is0 , T ie0l W is0 , 

^ F 3-4 frictional pressure loss across heat exchanger for 
isothermal conditions specified by T Jan 

and V it0 , lb/ft 2 ‘ 80 ' 

df 4 -5 frictional pressure loss through all ducts after 

heat exchanger up to final discharge section for 
isothermal conditions specified by „ 

ana w i g o • If desired, the pressure drop £P 4 _ 5 
may be subdivided into any number of smaller com- 
ponents, each of which must be corrected to non- 
isothermal conditions by the methods outlined in 
equation (6) of reference 78, lb/sq ft 

aF 5- 6 frictional pressure loss in isothermal discharge 
section for isothermal conditions specified by 
F isoj ^iso> an< * '^i so* 'lb/sq ft 

In equation (55) the terms on the left of the equal 
sign represent the difference in total pressure between the 
f.ee air stream and the point of air discharge. The first 
term on the right of the eq ual sign represents the pressure 

can easily be measured at sea level by an isothermal 
pressure loss test on a mock-up of the air distribution system. 


F is o 
^ F i - 2 

A F 3-3 
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changes due to the acceleration of the air in the duct, ow- 
ing both to changes in area and changes in specific volume. 
This term usually is quite small compared with the second. 
The second term represents the irrecoverable pressure loss 
due to the friction in the complete duct system. It should 
be noted that each isothermal frictional loss is corrected 
to- the operating temperature by different temperature cor- 
rections, depending or, the type of flow system represented 
by each separate AF. Thus, any complex flow system can be 
broken up into a series of systems, and the pressure drop 
through each corrected to noni sothermal conditions by the 
method outlined. The last term, APf an represents the to- 
tal pressure rise across a fan which may be placed in (say) 
the inlet duct to augment the ram pressure, The pressure 
change due to the adiabatic compression of the air between 
the free stream and the scoop entrance, sections 1-2, is 
neglected in this equation. This pressure change is small 
for usual aircraft speeds, but the temperature rise may be 
appreciable for velocities in excess of 3 00 ai.iles per hour 
and may be calculated from the equation ireference 78): 


where : 



k exponent for adiabatic compression in equation 
Pi v a k = V a k 

R gas constant for air, 53,3 ft-lb/lb °R 

In equation (55), for a given duct system for which the 
isothermal friction pressure loss AFj^g, AF 2 _ 3 , AF 3 _ 4 , AF 4 „ 5 

and AF 5 _ 6 are known, the remaining unknowns are W and T 4 

The fixing of the altitude, the airplane speed, and the heat 
loss from the duct establishes all other variables in the 
equation. Thus, for any altitude and airplane speed a curve 
of W versus T 4 can be drawn, which will reveal the rate 
of flow possible through the duct system for any temperature 

T 4 .- ' 


The relative importance of the various portions of the 
duct system may he readily established, for the largest of 
the corrected pressure drop terms in equation (55), will be 
the term which controls the rate of air flow. If it becomes 
necessary to increase the rate of flow through the heater- 
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duct system, attention should be focused on the largest term. 
By breaking up a complex duct system into a series of small 
units, the units causing a difficulty then may be readily 
isolated. 

Having established the curve of W versus T 4 from a 
consideration of the pressure drop characteristics of the 
duct system (from equation (55)), the thermal performance of 
the heater must be utilized in order to establish the operate 
ing point of the heater-duct system. The thermal performance 
of the heater is used to establish a second curve of W 
versus T 4 , which is fixed by the thermal output of the 
heater, since for any particular W and exhaust gas temper- 
ature, only one magnitude of T 4 is possible. The relation, 

(56) 


(57) 


is utilized to obtain this second curve. The heater capacity 
<llab' determined in the laboratory, must be corrected to 

altitude and temperature conditions by the method outlined in 
part III, section B-l, Temperature T 3 must include the 
temperature increase of the air due to compressibility be- 
tween, the free air stream and the scoop. The intersection 
of the curve of W versus T 4 obtained from the pressure 
drop characteristics of the heater— duct system (equation 
(55)), and the curve of W versus T 4 from equation (57), 
fixes the operating point of the system at the particular 
altitude and airplane speed under consideration, and allows 
the complete prediction of ventilating air rate, air temper- 
ature leaving the heater and heater output as a function of 
airplane speed and altitude. * 

For convenience in calculating the approximate perform- 
ance of the heater and duct system at various altitudes and 
airplane speeds, when flight data have Iseen obtained at one 
altitude and airplane speed the following simplifications of 
equation (55) are presented? 

*See reference J8 for a detailed example of such a pre- 
diction. 
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1. The first term on the right of the equal sign, which 
represents acceleration pressure drop, is usually 
negligible in comparison with the other terms. 
Thus, as a very close approximation (omitting the 
fan and the corresponding APf an ), 



2. If the velocity u 6 is very small and the static 
pressure P x is nearly equal to P 6 , which is 

the usual case for cabin heating, the equation 
reduces to: ¥ 


2 



2g ?! 


The 



ratio 


u- 


2g 


i s 


y 


a-h 



1 3 


(59) 


proportional to the square 


of the indicated airspeed Uj_ 2 . If, as a rough 

approximation, the temperatures of the air pass- 
ing through the duct are considered invariant with 
operating conditions, then the summation term is a 
constant for one duct-heater c omhination, Thus 


W 


k 




2 _ 

k(^ i/1?) * 


( 60 ) 


where k = constant, If the variation of the ex- 
ponent n from the second power is neglected, 
the rate of ventilating air flow due to the ram 

*If the velocity u 6 cannot he neglected, as may he the 
. . u ~ 2 

case m a wing anti-icirg system, the term must he 


retained in the equation 
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(when the velocity u 6 = 0) is approximately propor- 
tional to the first power o"f the indicated airspeed 
and the square root of the altitude pressure (when 
the velocity u 6 is small). If t>y test the rate of 
air flow through the heater and duct system is known 
at one altitude and indicated airspeed, the ventilat- 
ing air rate at any other altitude and airspeed can 
"be readily estimated by means of equation (60)*. Once 
the rate of air flow at various altitudes has been 
estimated, the temperature of the air leaving the 
exchanger may be calculated by means of the equation 


*alt " W a c p a ( T 4 - T s) 


(61) 


where 


q a it heater output at weight rate V , corrected 
to altitude conditions, Btu/hr 

V a ventilating-air rate, lb/hr 

Cp^ heat capacity of air, Btu/lb °F 

T 3 absolute temperature of air entering exchanger, 


absolxite temperature of air leaving exchanger, 




This method may be used to.obtain an estimate of 
heater and duct performance. The method is approxi- 
mate and reveals nothing of the pressure distribu- 
tion along the duct system. For a more detailed and 
precise analysis which will reveal the presstire dis- 
tribution, the graphical method presented in refer- 
ence J8 must be utilized. 

Referring again to equation (55), the following 
pertinent facts are noted. The terms of the form 



(62) 


represent the irrecoverable friction pressure loss 
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■between two sections of the flow system, a and h , 
for nonisothermal flow. It should he noted par- 


ticularly that tne temperature multiplier 


T a + T b y 
. 2T iso/ 


is greatest for the air just leaving the heater, be- 
cause, during normal operation, the air temperature 
is highest when just leaving the heater. Thus, 
special care must he taken to design the heater dis- 
charge section (or elhow) so as to make the isother- 
mal friction pressure loss £F a -b across this sec- 


tion as small £.s possible. If the isothermal fric- 
tion pressure loss across this section is large, the 
nonisothermalL friction pressure loss (due to the 
temperature multiplier) becomes excessive and may 
readily invalidate the advantages of a heater with 
low pres sure 'dr op. 


University of California, 

Berkeley, Calif., January 1944. 
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PROPERTIES OF AIR* 


T einnerature 


\x' x 10 ' 


(°F) (Btu/rb °F) (it sec/ft 3 ) ^Btu/hr ft 8 (^0)i 


Pi'andtl number 

= LfE (3600 g) 


-100 

0 

100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 


0. 2393 
.2398 
.2403 
.2412 
. 2427 
. 2449 
. 2476 
. 25 05 
. 2534 
.2566 
.2598 
.2630 
.2660 
.2690 
.2715 
.2740 
.2766 
.2789 


23 0 
343 
398 ' 5< ® 
449 * s 
498 
542 
587 '** 
63 0 ?0 J 
663^' \ 
699 3? ' 
732 ° 

7 67 H* 

8 00 ^ 
832 
864 
896 
928 
960 


0. 0104 
. 0130 
. 0157 
. 0182 
.0205 
. 0228 
. 025 0 
. 0272 
. 0293 
. 0314 
. 0334 
. 0355 
. 0376 
. 0399 
. 0419 
. 0440 
. 0461 
. 0484 


0. 743 -ft* 
. 731 
.706 
. 69 0 > C|1 
. 682 AW 
. 677 
. 672 
. 668 
. 666 
,6 63 
. 660 
. 658 
. 655 
. 652 
. 650 
. 648 
. 645 
.643 


*See references 17 and 79 for properties of other eases, 

la, Os, CO, COa, and H 3 . 







1 * 1 *] 
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A It i tude 
(ft) 

Temperature 
(°F ) 

Pressure 
(in. Hg.) 

Density 

(lb/ft 3 ) 

Density 

ratio 

0 

' 59.00 

29.92 

0.07651 

1.0000 

1,000 

55.43 

28.86 

.07430 

.9710 

2,000 

51.87 

27.82 

.07213 

.9428 

3,000 

48.30 

26.81 

.07001 

.9151 

4,000 

44.74 

25.84 

.06794 

.8881 

5,000 

41.17 

24.89 

.06592 

.8616 

6,000 

37.60 

23.98 

.06395 

.8358 

7,000 

34.04 

23.09 

.06202 

.8016 

8,000 

30.47 

22.22 

.06013 

.7859 

9,000 

26.90 

21.38 

.05829 

.7619 

10,000 

23.34 

20.58 

.05649 

.7384 

11,000- 

19.77 

19.79 

.05474 

.7154 

12,000 

16.21 

19.03 

.05303 

.6931 

13,000 

12.64 

18.29 

.05136 

.6712 

14,000 

9.07 

17.57 

.04973 

.6499 

15,000 

5.51 

16,88 

.04814 

.6291 

16,000 

1.94 

16.21 

.04658 

.6088 

17,000 

-1.63 

15.56 

.04507 

.5891 

18,000 

-5.19 

14.94 

.04359 

.5698 

19,000 

-8.76 

14.33 

.04216 

.5509 

20,000 

-12.32 

13.75 

.04075 

.5327 

21,000 * 

-15.89 

13.18 

.03938 

.5148 

22,000 

-19.46 

12.63 

.03806 

.4974 

23,000 

-23.02 

12.10* 

.03676 

.4805 

24,000 

-26.59 

11.59 

.03550 

.4640 

25,000 

-30.15 

11.10 

.03427 

.4480 

26,000 

-33.72 

10.62 

.03308 

.4523 

27,000 

-37.29 

10.16 

.03192 

.4171 

28,000 

-40.85 

9.72 

.03078 

.4023 

29,000 

-44.42 

9.29 

.02968 

.3869 

30,000 

-47.99 

8.68 

.02861 

.3740 

31,000 

-51.55 

8.40 

.02757 

.3603 

32,000 

-*55.12 

8.10 

.02656 

.3472 

33,000 

--58.68 

7.73 

.02558 

.3343 

34 ,000 

-62.25 

7.38 

.02463 

.3218 

35,000 

-65.82 

7.04 

.02369 

.3098 

30,000 

-67.00 

6.71 

.02265 

.2962 

37,000 

-67.00 

6.39 

.02 160 

,2824 

38,000 

-67.00 

6.10 

.02059 

• 2692 

39,000 

-67.00 

5.81 

.01963 

.2566 

40,000 

-67.00 ' 

5.54 

.01872 

.2447 

41,000 

-67.00 

5.28 

* .01785 

.2332 

42,000 

-67.00 

5.04 

.01701 

.2224 

43,000 

-67.00 

4.80 

.01622 

.2120 

44,000 

-67.00 

4.58 

.01546 

.2021 

45,000 

-67.00 

4.36 

.01474 

.1926 

46,000 

-67.00 

4.16 

.01405 

.1837 

47,000 

-67.00 

3.97 

.01339 

.1751 

48,000 

-67.00 

3.78 

.01277 

.1669 

49,000 

-67.00 

3.60 

.01217 

.1591 

50,000 

-67.00 

3.44 

.01161 

.1517 


*See reference 80 
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SUMMARY Of IQUATI0H8 


General form 


equations for air 


flat plate 


Laminar boundary layer 

, $ ° f x 

5 fPr} - £-• 


0.332 

3600 Uod y Op 3 V Re* 

Turbulent boundary l^yer 


3600 tt^r Op (Pr) “ 3 


Cf x 0.0396 


Re, 0 * 3 


0.5 

T 

0.5 

) 

Ol.5 

V x j 

1 

0.5 


\ 


V l , 



f e - 0.113 T f 
°av 1 

f c - 0.51 I f ° -3 (U " ^ 

c x 1 - 0.2 


0.8 


0.3 (Uoot) 
0*64 Tf TqJ 


0.8 


Pipes and duots 


Turbulent entranoe section 0 c x< 4.4 Dg 
*o T . 0.0396 

r^<"> - 7 ^ 

Beyond entranoe section; turbulent flow 
4.4 D H < x < oo 

a. 

3 t 0.176 


0 op 


“ Q " 


8 8 * Raj) 


0.3 


Laminar flow-parabolic velocity distri- 
bution at entrance 

(a) Round tubes 

f n 


- 1.16 ~\j 31 + 

k v kx 

(b) Root angular ducts 

f e ,<5 ?t 

-Sj- - 0.98 Y 59 + 


W Op 6 

kx B 


a 0.3 0°*B 
®x “ 73 * 10 ** ^3 


,_4 . 0.3 qO-8 
1.3 

*cav " ®- 1 x 10-4 T f°' 3 ?0 

X 


a 0.3 

f 0x - 5.4 X 10- 4 I ^ 


w av 


5.4 x 10' 


© 


4 0.3 Q 0,8 . . ®h\ 

T © 


-B 1 


f Cx » 3.65 \J\ ♦ 0.38 "y 

t 0% - 3.80 |\/l * 0.30 


flow aoross single cylinders 


1. Average over cylinder 

0.8 0,3 

HU - 0.26 Re Pr 


3. Local value at stagnation point (SP»0 ) 

„ , « 0-5 _0.4 

Hu • 1.14 Re Pr 

3. Local value along front half of 
cylinder (0<9 t 90°) 


Hu - 1.14 Re 


[> - a s ] 


1. 


3. 


3. 


f Q - 0.311 T 
°av 


0.« ( Uao y) 


. 0.6 


n0.4 


' w - 0 - 1941 ^' 49 ^) 

[ 


f^ - 0.194 T 


1 - -2- 
90° 


flow aoross tube banks 


Hu » 0 


0.6 | 

.244 F a Re Pr° 


o °** 

fo av - 14.5 X 10- 4 T f 0 - 43 Tt-frj 


Values of f- 


Humber of tubes 


3 


6 


8 


TO" 


In-line 
8t agger ed 


1.00 

1.00 


1.10 

1.11 


1.17 

1.33 


1.34 

1.31 


1.39 

1.39 


1.34 

1.45 


1.37 

1.48 


1.40 

1.51 


1.43 

1.53 


1.43 

1.54 


8t 

for Re ■ 20,000 and 1.35 < < 3.0 

1.35 < ^ < 3.0 
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